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stiffener  about  Its  line  of  attachment  to  the  plating  and  has  been 
demonstrated  by  some  recent  British  grillage  tests  to  have  serious 
potential  as  a  primary  mode  of  failure.  The  solutions  developed 
take  into  account  the  effects  of  (1)  the  rotational  resistance 
provided  by  the  plating  to  which  the  stiffener  is  attached,  (2)  non¬ 
linear  material  and  structural  behavior  by  means  of  a  tangent 
modulus  type  approach,  and  (3)  stiffener  web  deformations  (for  in¬ 
plane  loading).  The  equations  are  suitable  for  manual  calculations 
but  particularly  powerful  applications  are  possible  when  they  are 
teamed  up  with  a  desk-top  type  mini-computer. 

A  number  of  comparisons  between  tripping  predictions  made  using 
these  equations  and  numerical  finite  element  results  in  general  show 
very  good  agreement.  For  the  case  of  lateral  loading  the  agreement 
is  less  consistently  acceptable,  primarily  because  of  the  more 
complicated  nature  of  the  tripping  mode  shapes  and  the  inability  at 
present  to  include  the  effects  of  web  deformations  for  this  loading. 
Comparisons  are  also  made  with  experimental  collapse  data  from  two 
of  the  British  grillages  which  failed  by  tripping  and  the  agreement 
is  quite  good.  Unfortunately,  the  mean  tripping  failure  stresses 
(the  basis  of  these  comparisons)  are  quite  sensitive  to  the  esti¬ 
mates  of  plating  effectiveness  and  thus  wide  variations  in  predicted 
mean  stresses  are  possible.  Thus,  this  limited  experimental  vali¬ 
dation  is  not  as  conclusive  as  might  be  hoped  for. 

The  development  of  these  design  equations  is  quite  timely  in 
view  of  the  current  interest  being  shown  in  the  use  of  bulbs  and 
flat  bars  as  stiffening  members  (for  reasons  of  economy).  Such 
members  are  inherently  weak  with  respect  to  tripping  and  their 
application  would  require  that  careful  attention  be  given  to 
tripping  behavior. 
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ABSTRACT 


A  series  of  design  oriented  equations  to  predict  the 
tripping  (lateral-torsional  instability)  of  stiffeners 
under  inplane  and  lateral  loads  has  been  developed  which 
will  allow  this  mode  of  failure  to  be  more  comprehensively 
addressed  in  the  early  stages  of  structural  design.  This 
type  of  failure  is  characterized  by  a  twisting  of  the 
stiffener  about  its  line  of  attachment  to  the  plating  and 
has  been  demonstrated  by  some  recent  British  grillage 
tests  to  have  serious  potential  as  a  primary  mode  of 
failure.  The  solutions  developed  take  into  account  the 
effects  of  (1)  the  rotational  resistance  provided  by  the 
plating  to  which  the  stiffener  is  attached,  (2)  non¬ 
linear  material  and  structural  behavior  by  means  of  a 
tangent  modulus  type  approach,  and  O )  stiffener  web  de¬ 
formations  (for  inplane  loading) .  The  equations  are 
suitable  for  manual  calculations  but  particularly  power¬ 
ful  applications  are  possible  when  they  are  teamed  up 
with  a  desk-top  type  mini-computer. 

A  number  of  comparisons  between  tripping  predictions 
made  using  these  equations  and  numerical  finite  element 
results  in  general  show  very  good  agreement.  For  the 
case  of  lateral  loading  the  agreement  is  less  consistently 
acceptable,  primarily  because  of  the  more  complicated 
nature  of  the  tripping  mode  shapes  and  the  inability  at 
present  to  include  the  effects  of  web  deformations  for  this 
loading.  Comparisons  are  also  made  with  experimental 
collapse  data  from  two  of  the  British  grillages  which 
failed  by  tripping  and  the  agreement  is  quite  good.  Un¬ 
fortunately,  the  mean  tripping  failure  stresses  (the 
basis  of  these  comparisons)  are  quite  sensitive  to  the 
estimates  of  plating  effectiveness  and  thus  wide 
variations  in  predicted  mean  stresses  are  possible.  Thus, 
this  limited  experimental  validation  is  not  as  conclusive 
as  might  be  hoped  for. 

The  development  of  these  design  equations  is  quite 
timely  in  view  of  the  current  interest  being  shown  in 
the  use  of  bulbs  and  flat  bars  as  stiffening  members  (for 
reasons  of  economy) .  Such  members  are  inherently  weak 
with  respect  to  tripping  and  their  application  would 
require  that  careful  attention  be  given  to  tripping 
behavior. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  performed  at  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center  under  the  sponsorship  of  the 


NAVSEA  6.2  Exploratory  Development  Program,  specifically  the  Surface  Ship 
Structures  Block,  Program  Element  62543N,  Task  Area  SF  43  422  593  and 
Work  Units  1730-593  and  1730-610. 

INTRODUCTION 

A  series  of  tests  on  ship-type  steel  grillages  conducted  at  the  Naval 

1* 

Construction  Research  Establishment  (NCRE)  and  recently  reported  has 
clearly  demonstrated  the  significance  of  lateral-torsional  instability 
(tripping)  as  a  primary  ductile  failure  mode  for  ship  structure.  The 
potential  for  such  failure  has  important  ramifications  with  regard  to 
structural  weight,  fabrication  cost,  and  structural  reliability.  Whereas 
tripping  brackets,  which  are  perhaps  the  most  common  measure  employed  to 
prevent  this  ailment,  are  advantageous  from  a  weight  standpoint,  they 
suffer  the  disadvantages  of  increasing  fabrication  cost  and  of  introducing 
hard  spots  at  their  toes  which  may  give  rise  to  fatigue  problems  and  weaken 
the  structure  under  explosive  loading.  In  any  case  it  seems  clear  that  the 
ability  to  predict  tripping  failure  in  the  early  stages  of  design  can  have 
important  consequences  on  the  design  of  both  conventional  and  high 
performance  ships. 

Surprisingly  little  material  exists  in  the  literature  on  the  subject 

of  tripping  of  stiffeners  welded  to  continuous  plating.  In  the  elastic 

2  3 

region  tripping  stresses  can  be  estimated  using  approximate  formulas  ’ 

(for  certain  selected  cases  of  loading)  or  more  generally  using  folded- 
4  5 

plate  *  or  finite  element  analysis.  In  the  inelastic  region  no  satis¬ 
factory  method  appears  to  exist  at  present  although  the  application  of 
incremental  finite  element  analysis  would  appear  to  offer  great  promise. 

The  main  thrust  of  this  study  is  the  development  of  fast,  approximate 
methods  of  tripping  analysis  for  particular  application  in  the  concept, 
feasibility,  and  preliminary  stages  of  design.  As  mentioned  above, 
approximate  formulas  do  exist  in  the  literature,  however,  the  work 
described  herein  both  modifies  and  extends  these  previous  efforts.  As 
with  these  previous  efforts,  the  solutions  developed  are  elastic  in  nature; 
inelastic  effects  must  be  accounted  for  by  modifying  the  elastic  solutions, 


*A  complete  listing  of  references  is  given  on  page  87. 


an  admittedly  crude  approach.  Nevertheless,  the  value  ol  having  re  I  at lvel \ 
simple,  approximate  solutions  at  the  early  design  stages  fwhcn  more  sophis¬ 
ticated  solutions  are  otten  not  practical)  cannot  he  overstated;  many  ol 
the  most  significant  and  critical  decisions  are  made  during  these  early 
design  stages. 

TRIPIMNi;  UNDER  AXIAI.  OK  IjVTK.RAl.  LOADS 
l.ateral-torslonal,  or  tripping  instability  is  characterized  hv  a  twist 
Ittg  of  the  stiffener  about  Its  line  ot  attachment  to  the  plating.  This 
deformation  pattern  Involves  both  sideways  and  vertical  flexure  (v,w)  and 
rotation  (t')  ot  the  stiffener  as  shown  in  Figure  I.  It  the  centerline  ot 


0 

Figure  la  -  Coordinate  System 


Figure  l  -  Character izat ion  of  Tripping 


the  web  at  any  transverse  cross  section  is  not  allowed  to  curve,  then  the 
three  displacement  degrees  of  freedom  v,  w,  and  3  are  coupled  since,  for 
small  displacements,  v  %  z6  and  w  «  -y6. 

The  critical  load  for  tripping  Is  approximated  by  applying  Rayleigh's 
principle.  This  principle  states  that  for  all  possible  deformations,  the 
total  strain  energy  V  is  greater  than  or  equal  to  the  work  W  done  by  the 
externally  applied  forces  (or  by  the  internal  stress  field  arising  from 
such  forces).  As  a  consequence,  it  is  possible  to  obtain  an  approximate 
critical  buckling  stress,  one  that  is  always  greater  than  the  true 
buckling  value,  by  equating  the  total  strain  energy  and  external  work, 
l.e. , 


V  -  W  (1) 

Inherent  in  the  method  is  the  assumption  of  a  buckled  mode  shape;  equating 
the  energies  as  indicated  by  Equation  (1)  will  provide  either  upper  bounds 
to  the  critical  stresses  or  exact  solutions  in  the  unlikely  event  of  the 
assumed  mode  shape  being  exact. 

The  total  strain  energy  stored  in  the  structure  in  the  buckled  state 
(neglecting  the  energy  of  vertical  bending  in  the  flange)  is  given  by  the 
following  expression 

a 

V  -  ~  I"  [El  v2  +  Er  £2  +  GJ  ft2  +  C62 1  dx  (2) 

2  J  z  xx  xx  x 

0 

where  v  and  6  here  refer  to  the  translation  of  and  rotation  about  the 
shear  center  of  the  stiffener  and  the  x  subscripts  indicate  partial 
derivatives  with  respect  to  x.  The  strain  energy  is  made  up  of  contri¬ 
butions  from  (in  the  order  of  appearance  above)  sideways  bending,  longi¬ 
tudinal  warping,  torsion,  and  rotation  of  the  supporting  plate  structure 
modeled  as  an  elastic  spring.  Since  the  sideways  flexure  of  the  shear 
center  is  coupled  to  the  rotation  P  as  indicated. 


4 


V  ~  80 


(3) 


the  sideways  bending  and  longitudinal  warping  terms  can  be  combined  to 
give 


a 

V  -  |  |  [E(IzS2+D  02x  +  GJ  02  +  C02]  dx  (A) 

0 

In  this  expression  a  Is  the  length  of  the  stiffener  between  transverse 
supports,  E  and  G  are  the  Young's  and  shear  moduli  of  the  material,  re¬ 
spectively,  and  C  Is  the  rotational  spring  constant  (per  unit  length)  of 
the  supporting  structure.  The  coefficients  1^,  s,  P,  and  J  are  defined  in 
Figure  2  in  terms  of  the  geometry  of  the  stiffener. 

Although  Equation  (A)  may  be  applied  to  flat  bars  as  well  as  flanged 
stiffeners,  it  is  perhaps  more  appropriate  for  thin  flat  bars  to  use  the 
expression  for  strain  energy  derived  from  thin  plate  theory,  namely 

a  d(x) 

V  *  t  D  [  f  {(v+v)2-  2(l-v)  [v  v  -  v2  ] }  dxdz 
2  w  J  J  xx  zz  xx  zz  xz 

0  0 

a 

+  \  j  C02  dx  (5) 

0 

where  is  the  well  known  plate  flexural  rigidity  (for  the  flat  bar), 

V  is  Poisson's  ratio,  and  d(x)  represents  the  depth  of  the  flat  bar  at 
location  x. 

In  computing  the  work  W  the  nature  of  the  external  loading  will  often 
dictate  the  form  of  the  expressions  used.  When  axial  end  thrusts  and/or 
moments  are  applied  at  the  stiffener  ends  it  is  usually  most  convenient  to 
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Figure  2 


lz  -  Movent  of  Inertia  about  the  Ueb  Plane 
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s  -  Height  of  Shear  Center  above  Toe  (Origin) 
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r  -  Longitudinal  Harping  Constant 
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It  -  Vertical  Moment  of  Inertia  about  Toe  (Stiffener  alone) 

I  -  i  t  d3  +  f  tf  (d2  +  7Y  t3) 
t  3ww  wfvc  12  f  / 


I  -  Polar  Movent  of  Inertia  about  Toe 
P 


z  -  Height  of  Centroid  above  Toe 

z  -  U  t  d2  +  f  t,d  1  /ft  d  +  f  t.l 
L 2  w  w  w  f  CJ/  [  w  v  w  f  J 

-  Geometrical  Tripping  Parameters  for  Tee  Stiffeners 
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compute  the  component  of  work  done  by  these  forces  by  integrating  the 
product  of  the  axial  stresses  arising  from  these  forces  and  the  axial 
shortening  6(y,z)  of  the  stiffener  at  the  point  in  question,  namely. 


" '  JJ°> 


6(y,z)  dydz 


(6) 


As  the  expression  suggests,  the  integration  is  carried  out  over  the  area 
of  the  stiffener  end.  When  the  stiffener  has  lateral  loading  the  above 
approach  cannot  be  applied,  rather  the  work  done  by  the  internal  stress 
field  arising  from  the  lateral  loading  is  computed  according  to 


W 


1 

2 


d 

a  1 

J  / 

0  0 


[a  (v  )2  +  o  (v  )2 
X  X  z  z 


2  T  v  v  ]  t  dzdx 
xz  x  z  w 


C 

J 


[o  (v  +w  ) ]  .  f  t  dx 

x  x  x  z-d  w  f 
c 


(7) 


where  the  first  integral  represents  the  action  of  the  stress  field  in  the 
web  and  the  second  that  in  the  flange.  When  both  lateral  and  end  loadings 
are  present,  the  stress  fields  arising  from  each  set  of  loads  may  be 
separated  and  the  corresponding  work  components  computed  according  to 
Equations  (6)  and  (7),  as  appropriate. 

The  signs  of  the  terms  in  Equation  (7)  differ  in  some  respects  from 
those  found  in  other  references.  This  is  due  to  the  assumption  of 
compression  as  a  positive  stress  and  tension  as  negative.  (This  is  quite 
common  in  compression  instability  work.)  Figure  3  defines  the  sign 
convention  adopted  in  this  report  for  the  o^,  o^,  and  Txz  stresses  and  the 
corresponding  equilibrium  equations  which  are  consistent  with  Equation  (7). 
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EQUILIBRIUM  EQUATIONS 
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do 

_ * 

9z 
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Figure  3  -  Sign  Convention  for  Inplane  Axial  and  Shear 
Stress  Components 


TRIPPING  UNDER  END  LOADS 

Consider  a  single  stiffener  of  constant  cross  section  and  its 
associated  frame  space  of  plating  loaded  by  an  axial  inplane  force  of 
magnitude  P  as  shown  in  Figure  4.  The  ends  of  the  stiffener  will  ex¬ 
perience  a  uniform  axial  stress  o^  which  is  related  to  the  cross  sectional 
area  of  the  stiffener  A^,  the  plating  thickness  t,  and  the  plating 
effective  width  be,  as  follows. 


o 

e 


P 

A  +  b  t 


s  e 


(8) 


In  terms  of  the  stress  components  in  the  stiffener,  o^,  o^,  and  T  ,  this 
translates  into 


Figure  4b  -  Constant  Moment 


T  I  I  I  t  I 


■ti 


Figure  4c  -  Uniform  Lateral  Load 
Figure  4  -  Geometry  of  External  Loadings 


a 
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a 
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(9) 


a 

z 


0 


(10) 


To  calculate  the  strain  energy  in  the  stiffener  and  the  work  done  by 
the  axial  load,  an  assumption  has  to  be  made  for  the  buckled  mode  shape. 
Because  of  the  repeatability  and  continuity  of  typical  ship  structure,  it 
is  logical  to  select  a  mode  shape  incorporating  simply  supported 
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boundaries,  (6^*0)  at  the  ends  x  ■  0  and  x  ■  a  In  addition  to  the  require¬ 
ments  for  zero  displacements  (6*0) .  Perhaps  the  simplest  choice  which 
satisfies  these  conditions  is  to  assume 
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■ 


which,  in  turn,  upon  substitution  of  Equation  (12)  becomes 

<5(y.z)  “  £  I*"  +  y2]  (;p)  a  <17> 

The  work,  therefore,  is 

W  “  °e  *  i  “  JJ"  lz"  +  y2]  dydz  (18) 

A 

s 

where  the  integration  over  the  stiffener  cross  section  is  carried  out  at 
the  point  of  load  application.  Since 

[z“  +  y~]  dydz  “  1  (19) 

s 


(see  Figure  2)  the  expression  for  the  work  of  the  external  force  may  be 
further  transformed. 


w-±o  a2  (51) 

A  e  ova' 


a  I 


(20) 


Equating  the  total  strain  energy  and  the  work  leads  to  the  critical 
elastic  value  of  stress  for  stiffener  tripping.  For  flanged  stiffeners 
this  stress  is 


(°e>cre "  (0 


and  for  flat  bars. 
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(o  ) 

e  ere 


D  L  J  D 


where  (for  flat  bars) 


t  .  1  ,  ,3  4.  1  .,.3  1  ,3 

I  ■  -r  t  d  +  at  ®  -=■  t  d 

p  3  w  12  w  3  w 


In  actual  fact  these  solutions  are  exact,  as  comparisons  with 

published  solutions**  will  demonstrate.  (This  is  of  course  subject  to  the 

condition  stated  earlier  of  zero  curvature  v  of  the  stiffener  web.) 

zz 

Thus,  Kquation  (12)  in  this  case  represents  the  exact  tripping  mode  shape. 

When  there  is  no  restraint  against  rotation,  C  -  0,  it  can  be  seen  by 

inspection  that  the  lowest  buckling  stress  occurs  for  one  wave,  m  ■  1. 

When  C  is  nonzero,  the  mode  number  for  which  (o  )  is  lowest  is 
’  e  ere 

dependent  on  the  degree  of  restraint.  Although  m  takes  on  only  integer 
values,  the  determination  of  the  critical  stress  can  be  simplified  by 
setting 


3(0  ) 
e  ere 


solving  for  m,  and  then  comparing  the  tripping  stresses  for  the  two  integer 
values  of  m  which  bracket  this  value.  From  the  above  conditions  one 
obtains 


.  *  r _ £ _ 1 

17  [_e(i  s2+n  J 


for  flanged  stiffeners,  and 


a  r  3C 
71  Dd3 


for  flat  bars. 
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e  ere 

at  the  stiffener  location  and  tl»at  when  the  effective  width  of  the  plating 

bg  is  less  than  b  this  stress  will  be  greater  than  both  the  corresponding 

average  stress  in  the  plating  o  and  the  mean  stress  over  the  total  plate- 

beam  cross  section,  0  .  These  definitions  of  stress  are  illustrated  in 

m 

Figure  5.  Since  the  mean  stress  is  frequently  used  in  design  calculations 


Figure  5  -  Distribution  of  Inplane  Stress 


the  peak  stresses  given  by  Equations  (21)  and  (22)  can  readily  be  con¬ 
verted  to  the  mean  form  according  to 

/  A  yfb  t  \ 

(0.>cr.-(in5r)<0.>cr.  (27> 

'  s 

The  distinction  between  peak,  average  (plate),  and  mean  stresses  is  a 

particularly  significant  point  with  regard  to  tripping  since  it  is  the 

stress  which  actually  occurs  in  the  stiffener  which  is  critical  and  the 

treatment  of  (°e)cre  as  either  an  average  or  a  mean  stress  can  lead  to 

very  unconservative  estimates  of  the  loads  required  to  cause  tripping. 

For  example,  current  Navy  design  practice  assumes  that  primary  stresses 

are  distributed  uniformly  over  the  total  plate-beam  cross  section  provided 

that  these  stresses  are  less  than  those  required  to  cause  plate  buckling. 

If  this  is  true  (and  only  primary  stresses  are  present)  then  all  the 

stresses  are  identical,  i.e.,  o  *  o  «  o  .  In  actual  fact,  however, 

earn 

imperfections  associated  with  real  structure  will  cause  a  reduction  in 

plating  effectiveness  even  prior  to  plate  buckling  and  lead  to  a  nonuniform 

distribution  of  primary  stress  as  Figure  5  illustrates.  In  this  situation 

the  failure  to  take  into  account  the  fact  that  the  stress  in  the  stiffener 

is  higher  than  Lhe  calculated  mean  stress  could  lead  to  an  unpleasant 

surprise  since  the  tripping  stress  in  the  stiffener  will  be  achieved  at  a 

loading  somewhat  lower  than  that  predicted  by  the  assumption  of  uniform 

primary  stress.  The  degree  to  which  peak  and  mean  stresses  will  differ  is 

a  function  of  plating  effectiveness,  as  represented  by  the  effective  width 

(and  the  effective  breadth  where  secondary  loads  are  also  present)  and  the 

accurate  prediction  of  stiffener  tripping  is  clearly  dependent  on  the 

ability  to  reasonably  estimate  this  effectiveness. 

If  the  plate-stiffener  combination  is  loaded  by  a  constant  moment  M, 

as  show  in  Figure  4b,  the  stress  distribution  will  no  longer  be  constant 

over  the  cross  section  but  rather  will  vary  linearly  from  plate  to  flange. 

Thus  in  terms  of  0  ,  0  ,  and  X  one  obtains 
X  z  xz 

°x  •  (f )  H*  1 0  <28> 
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0 


(29) 


T  -  0  (30) 

xz 

where  I  is  the  vertical  moment  of  inertia  of  the  stiffener  and  its 
associated  effective  breadth  of  plating  and  h  is  the  corresponding  location 
of  the  neutral  axis  from  the  midplane  of  the  plating. 

If  the  same  mode  shape.  Equation  (12),  is  assumed,  then  the  ex¬ 
pressions  for  total  strain  energy  are  once  again  given  by  Equations  (13) 
and  (14).  The  work  is  computed  using  Equations  (6)  and  (17)  leading  to 


“-£eo(?)2  ■(?)// 


o  2 

[*“  +  y  1 


h  +  -7  t 


dydz  (31) 


The  integral  in  this  expression  is  broken  up  into  two  parts,  one  for  the 
web  and  one  for  the  flange 


if  [z2  +  y2)  [2  -  h  +  y  tj  dydz 
* JA 


r  w 

.f 

2  w 

[ 

2  2 
[z  +  y  ] 

^z  -  h  +  ~  t  j 

0 

-  A  t 

2  w 

r  j 

f  i  f 
.  2  w 

.2  2 . 

[z  +  y  ] 

z  -  h  +  t  j  ( 

d 

w  - 

\  f 

2  w 

dydz  (32) 
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Because  the  web  depth  and  flange  width  are  typically  much  larger  than  the 
web  and  flange  thicknesses,  these  integrals  can  be  simplified  by  making 
certain  approximations,  resulting  in 


fj  [z2  +  y2]  -  h  +  j  tj  dydz 

A 

s 

*  “  (h  '  \  C)  JJ  (z2+y2)  dyds 


I 


+  K 


*f  |  z  t  dz  + 
w 


I 


dc(dc+y2)  cfdy 


-  1  fW 


-  (h-  |  t)  I  +  \  t  dA  +  d  (d2t  f  +  yy  f3t  ) 
\  2/p  4ww  c\cfw  12  w  f  / 


td4-4(h-yt')l  +  4  d  (Afd2+I 

ww  \  2/p  c  fczf 


(33) 


where  I  and  A^  are  the  moment  of  inertia  about  the  web  plane  and  the 
cross  sectional  area  of  the  stiffener  flange  alone,  respectively.  Because 
this  geometrical  parameter  will  continually  reappear  in  the  later  sections 
of  this  report  it  has  been  denoted  by  the  symbol  S  as  indicated  above. 

The  above  expression  is  valid  for  both  flanged  stiffeners  and  flat  bars 
although  for  flat  bars  the  last  term  is  obviously  zero  since  both  A^  and 
I  ^  are  identically  zero. 

Making  use  of  this  definition,  the  work  done  by  the  applied  moment 
is  given  by 

"-f  eo(;r)2a(i)s  (34) 
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As  with  the  case  of  a  concentrated  force,  the  critical  moment  for 
elastic  tripping  is  determined  by  equating  the  total  strain  energy  and 
the  work.  The  resulting  expressions  are 


M 


ere 


E<V2+r>  (?) 


2 


+  GJ  + 


2 


C 


(35) 


for  flanged  stiffeners,  and 


M 


ere 


4  D  d 
3  w 


{(?) 


d2  + 


6(l-v) 


(-) 

\m  / 


(36) 


for  flat  bars. 

For  the  case  of  no  restraint  against  rotation  (C=0)  the  critical 
moment  also  occurs  for  the  mode  shape  corresponding  to  m  =  1.  When  re¬ 
straint-  does  exist,  the  critical  moment  may  be  found  using  the  same  pro¬ 
cedure  as  that  described  for  constant  force  loading,  and  in  fact,  it  is 
quite  easy  to  show  that  the  values  of  m,  for  which  Mcrg  is  minimum,  are 
given  by  the  identical  expressions,  Equations  (25)  and  (26).  It  is  also 
noteworthy  that,  for  all  practical  considerations,  only  positive  moments 
will  cause  tripping  as  indicated  by  Equations  (35)  and  (36).  Only  when  the 
neutral  axis  is  located  at  a  significant  distance  from  the  plating  may  a 
negative  moment  be  likely  to  cause  tripping  and,  in  this  case,  the  assumed 
deformation  used  in  the  preceding  development  would  probably  not  be 
appropriate. 


TRIPPING  UNDER  UNIFORM  LATERAL  LOADING 

Many  stiffened  plate  components  of  a  ship  routinely  experience  lateral 
loadings  such  that  the  possibility  of  tripping  collapse  due  to  such  load¬ 
ings  (or  at  least  in  conjunction  with  inplane  loading)  is  a  definite 
possibility.  If  this  type  of  failure  is  to  be  avoided,  it  is  clear  that 
the  ability  to  predict  tripping  behavior  under  lateral  loading  is 


17 


necessary.  Since  in  many  design  calculations  lateral  loadings  are 
idealized  as  uniformly  distributed,  the  development  which  follows  assumes 
such  a  distribution.  However,  the  technique  employed  is  not  restricted  to 
uniform  distributions  and  can  readily  be  applied  to  others  as  the  need 
arises . 

The  critical  lateral  load  for  stiffener  tripping  is  computed  in  the 
same  manner  as  that  for  concentrated  end  loading,  that  is,  the  strain 
energy  V  is  equated  with  the  external  work  W  for  an  assumed  displacement 
function.  In  this  case,  however,  W  must  be  determined  by  computing  the 
work  done  by  the  internal  stress  field  arising  from  the  presence  of  the 
lateral  loading  as  indicated  in  Equation  (7) . 

For  a  single  stiffener  and  its  associated  frame  spacing  of  plating, 
simply  supported  at  its  ends,  the  stress  components  due  to  a  uniform 
lateral  loading  q,  (positive  in  the  positive  z  direction)  are  as  follows. 


°X  =  i(i)  ( z~h+i c)  (x2-ax) 

/  \  /  rA,z  .  .  . 

(f  -(f)  r-  k--T*) 

\  w  /  L  w 


+  2  {i  dv  -  i  "2  +  (h~  1 2)  (i2-d„)| 


\z-i(?)(a-2lt>fr(dc-h+i 


}d2-}z2  +  U  t )  rz-d  )1 

2  w  2  \  2  /  w 


(These  stresses  can  readily  be  shown  to  satisfy  the  equilibrium  equations 
given  in  Figure  3.) 


Although  «WU  appear  logical  to  again  adopt  lauation  1Z>  as  tne 
assumed  displacement  function,  -experience  has  Indicated  char.  lue  to  cue 
concentration  of  high  CT  stresses  around  the  uldlangth  of  -he  stiffener 
this  function  is  in  serious  error  with  re  sneer  to  the  true  tic  idea  itate. 
An  improved  displacement  function  tan  he  onstructed  by  nr  I noi-ig  i  compo¬ 
nent  of  the  next  higher  mode  along  with  the  primary  node,  sneer f -Trail vT 

If  sin  r~  r  (K— 1)  sin  tStQJZ  ,40) 


where  m  is  the  primary  mode  number.  This  expression  incorporates  an 
unknown  coefficient  if,  which,  in  general,  nay  take  on  values  oetveen 
approximately  1.5  and  0.  Since  the  application  of  Sayleigh 1 3  principle 
produces  a  htickling  load  which  is  always  greater  than  the  true  value,  the 
appropriate  value  of  if  is  selected  by  minimizing  the  value  of  trie  buckling 
load  with  respect  to  if. 

The  strain  energy  expressions  are  determined  by  substituting  equation 
(40)  into  Equations  (4)  and  (5)  and  performing  the  necessary  operations. 
This  leads  to 


r'z,4(z)  tU) 

for  flanged  stiffeners,  and 

2 

V  -  ~  t>  ad  p,2  (  -  )  H  (K)  (42) 

12  w  o'a'm 


for  flat  hats.  In  these  expressions  II  (K)  and  H  (K)  are  shorthand 

m  m 

notation  fot  quadratic  functions  in  K  which  are  defined  as  follows: 
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Hm(K)  =  (l1)2  ^  +  (™+2)4  (K-l)2} 


+  6(l-v)  {m2K2  +  (m+2) 2  (K-l)2} 

+  Fd(f)  {K2  +  (K-l)2} 
w 

The  work  W  done  by  the  internal  stress  field  is  calculated  using 
Equation  (7),  giving 

W  ‘  '  *8  (?)  6o  "2  *  S  F„<«  («> 

where  once  again,  for  reasons  of  convenience,  shorthand  notation  is  used 
to  represent  a  quadratic  function  in  K,  namely 

V«  '  "2 1(2  f1-  -tt)  +  *  1 1 

\  m  TT  /  TT  y 


+  ■  V  K(l-K)  +  (m+2)2  (K-l)2  1 

(m+1) ^  l 


Equating  the  strain  energy  and  the  work,  and  then  solving  for  q 
results  in  the  critical  value  of  the  lateral  loading.  For  flanged 
stiffeners  one  gets 


-  _  12IGJ  #  Hm(K) 

qcre  “  2  *  F  (K)  (46> 

ba  in 


and  in  a  similar  fashion  for  flat  bars. 


4ldDw  H  <K) 

qcre  =  2  '  F  (K)  <47) 

Sa  m' 


The  above  expressions  for  qcrg  contain  the  parameter  K.  whose  value  at 

this  point  is  still  unknown.  As  previously  mentioned  the  appropriate 

value  of  K  is  determined  by  minimizing  the  value  of  q  with  respect  to  K. 

ere 

This  process  is  routinely  carried  out  by  setting  the  derivative, 

3qcre/3K  =  0  and  solving  for  K. 

In  the  above  expressions  for  qcre,  the  K  terms  only  appear  in  the 

and  H^lO/F^K)  ratios.  Since  the  H's  and  F's  are  quadratic 
functions,  they  can  symbolically  be  represented  in  the  form 

VK)  =  h0  +  hlK  +  H2r2  (48a) 

Fm(K)  “  f0  +  flK  +  f 2r2  (48b) 

It  can  readily  be  shown  that  setting  the  derivative  9q  /3K  equal  to  zero 

ere 

leads  to  a  third  quadratic  equation  in  K, 

lhlf2  ~  Vl*2  +  2[h0f2  -  h2f0]K  +  [h0fl  "  hlf0]  *  0  (49) 


21 


the  solution  of  which  provides  the  value  of  K  for  which  H  (K)/F  (K),  and 

m  m 

hence,  <lcre>  is  a  minimum.  The  above  equation  being  quadratic,  usually 

two  distinct  values  of  K  will  result.  The  appropriate  choice  of  K, 

however,  is  usually  quite  obvious  as  is  discussed  later  in  Appendix  A  in 

connection  with  the  illustrative  examples. 

The  expressions  for  H  (K) ,  H  (K) ,  and  F  (K)  as  given  by  Equations  (43) 

mm  m 

and  (45)  are  not  in  the  form  required  for  the  determination  of  K  according 
to  the  method  represented  by  Equation  (49).  These  expressions  have  been 
rearranged,  however,  and  the  forms  most  appropriate  for  use  with  Equation 
(49)  are  presented  in  Appendix  B. 

If  the  ends  of  the  stiffener  and  its  associated  frame  space  of  plating 
are  clamped  rather  than  simply-supported,  the  stress  components  due  to  a 
uniform  lateral  loading  q  are  as  follows, 

°x  ■  7  (  I )  <1‘Z-ax)  (z‘h+  1  0  +  II  (?)  “2  (z-h+  1  ')  <50) 


+  z 


{K-S-’  +  M*)  (i«0) 


(51) 


xz  2 


I  /V 


(!) (a_2x)  IT  (dc"h+  1  c) 


V-  W 


w 


1  2  . 
2Z  + 


(z-d  ) 
w 


(52) 


If  Equation  (40)  is  again  selected  as  the  assumed  displacement 
function,  then  the  strain  energy  expressions  are  the  same  as  for  the 
simply  supported  case,  namely  Equations  (41)  and  (42).  The  expressions 
for  the  external  work,  however,  will  not  be  Identical  since  the  stress 
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component  is  different.  The  differences  are  entirely  within  the  definition 

of  the  function  F  (K),  however,  so  that  Equation  (44)  is  valid  for  the 
m 

clamped  condition  with  F  (K)  now  defined  as 

m 


F  (K) 
m 


K  K2  +  ~  (m+2)  Jl  +  -1— -=  K(l-K)  -  (K-l) 2 

n  7T  l  (m+l)^)  TT 


6d 

w 

7T2S 


Aftf 


(dc-M  I  C) 


+  b  t 
e 


(K2  + 


(K-l)2] 


(53) 


Using  this  definition  of  Fm(K)  the  critical  value  qcre  may  be  computed 
using  Equations  (46)  and  (47)  and  the  procedure  previously  outlined  to 
determine  the  appropriate  value  of  the  parameter  K.  The  expressions  for 
the  "clamped"  version  of  Fm(K)  in  the  format  which  is  symbolically 
represented  by  Equations  (48)  is  also  provided  in  Appendix  B. 

COMPARISONS  WITH  NUMERICAL  SOLUTIONS 

To  demonstrate  the  application  of  these  solutions  and  to  shed  some 
light  on  their  approximate  ranges  of  validity,  several  comparisons  were 
made  with  solutions  obtained  by  the  finite  element  method,  specifically 
with  a  program  developed  and  documented  at  the  University  of  California. ^ 
Both  a  flat  bar  and  a  tee  stiffener  configuration  were  considered,  the 
proportions  of  which  are  provided  in  Figure  6. 

Because  of  the  symmetry  of  the  structure  and  its  response,  it  was 
necessary  to  model  only  one-half  of  each  stiffener  for  the  finite  element 
analyses.  A  uniform  mesh  consisting  of  10  rectangular  plate  elements  in 
the  x-direction  and  8  elements  in  the  z-direction  was  employed.  The  flange 
of  the  tee  stiffener  was  modeled  by  an  additional  10  bar-type  elements. 
Although  assumptions  inherent  to  the  finite  element  method  mean  that  so¬ 
lutions  obtained  using  the  method  are  themselves  approximate,  the  moder¬ 
ately  fine  mesh  employed  here  can  be  expected  to  produce  results  which  are 
at  most  a  few  percent  different  from  the  so-called  theoretically  "exact" 
values. 
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TEE  STIFFENER  -  160  cm 

Figure  6  -  Stiffener  Geometries  for  Comparative  Solutions 

Tables  1  and  2  provide  buckling  stress,  moment,  and  lateral  load 
coefficients  for  the  first  two  modes  for  the  case  of  zero  rotational  re¬ 
straint  at  the  stiffener  base,  i.e.,  C  =  0.  For  the  cases  of  axial  end 
load  and  end  moment  agreement  is,  in  general,  very  good  between  the  finite 
element  and  energy  solutions.  For  lateral  loading,  the  agreement  between 
the  two  methods  of  solution  varies  from  very  good  to  very  poor,  with  the 
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TABLE  1  -  TRIPPING  STRESS  COEFFICIENTS  UNDER  END  LOADS 
AND  MOMENTS  (NO  ROTATIONAL  RESTRAINT) 


Flat  Bar 

Mode 

m 

Tripping  Stress  [(o  )  /E]  x  10^ 

Coefficient 

Finite  Element 

Eqn.  (22) 

i 

6.44 

6.42 

2 

7.58 

7.51 

Mode 

m 

Tripping  Moment  ("m  (d-h+  ^  t)  x  10^1  /  [IE] 

Coefficient  L  cre  '  '  J  ' 

Finite  Element 

Eqn.  (36) 

i 

9.04 

8.99 

2 

10.66 

10.52 

Tee  Stiffener 

Mode 

m 

Tripping  Stress  I(oe)cre/E]  x 

Coefficient 

Finite  Element 

Eqn.  (21) 

i 

26.4 

26.3 

2 

65.9 

72.6 

Mode 

m 

Tripping  Moment  M  (d  -h+  y  t)  x  10^  /  [IE] 

Coefficient  L  Cre  c  i  /  J 

Finite  Element 

Eqn.  (35) 

i 

29.4 

29.0 

2 

75.1 

80.1 
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TABLE  2  -  TRIPPING  LOAD  COEFFICIENTS  UNDER  UNIFORM 
LATERAL  LOAD  (NO  ROTATIONAL  RESTRAINT) 


Flat  Bar 

Mode 

m 

B.C.* 

a2(d-h+  j  t) 

Tripping  Load  100  q 

Coefficient 

Finite  Element 

Eqn.  (47) 

1 

2 

SS 

SS 

-  1.13 

-  1.20 

-  1.12 

-  1.18 

1 

2 

C 

C 

+  1.76 

+  3.18 

+  1.96 

+  3.90 

Tee  Stiffener 

Mode 

m 

B.C.* 

a2(dc-h+  i  t) 

Tripping  Load  100  q 

Coefficient  1 

Finite  Element 

Eqn.  (46) 

1 

2 

SS 

SS 

-  4.59 

-  9.03 

-  5.20 

-10.09 

1 

2 

C 

C 

+  6.54 

+17.91 

+  7.91 

+40.76 

*B.C.  =  Boundary  Conditions,  SS  =  Simple 

Support,  C  =  Clamped. 
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flat  bar  simply  supported  showing  the  best  agreement  and  the  clamped  tee 
stiffener  the  worst.  If  restraint  is  added  along  the  base  of  the  stiffen¬ 
ers  the  discrepancies  between  the  two  methods  consistently  grow  as  the 
amount  of  the  restraint  is  increased.  This  trend  is  clearly  illustrated 
in  Tables  3,  4,  and  5  where  buckling  coefficients  are  tabulated  as  a 

function  of  a  dimensionless  restraint  coefficient  R  ,  defined  as 

c 


R 

c 


gj  +  E<i2s2+n  (si) 


2 


(54) 


for  flanged  stiffeners,  and 


R 

c 


4  D  d 
3  w 


<k)l 

(j1)  d2  +  6(l-v) 


(55) 


for  flat  bars. 

If  the  finite  element  data  is  examined  in  more  detail  it  becomes 

quite  apparent  that  the  primary  cause  of  the  discrepancies  noted  above  is 

the  presence  of  web  deformations.  To  prevent  overly  optimistic  buckling 

3 

predictions,  Faulkner  in  his  work  suggested  an  upper  limit  on  the  value 
of  C  used  as  follows, 


Ca 


2 

ir  GJ 


<  10 


(56) 


For  the  tee  stiffener  cited  here,  however,  this  corresponds  approximately 
to  a  value  of  Rc  =  4.1  (for  m=l)  and  from  the  data  presented  in  the 
tables,  the  above  limit  itself  would  appear  to  be  overly  optimistic. 

While  the  lack  of  consideration  of  web  deformations  appears  to  be 
the  primary  cause  of  the  discrepancies  noted  in  the  tables,  a  possibly 


27 


TABLE  3  -  TRIPPING  STRESS  COEFFICIENTS  FOR  FLAT  BARS  UNDER  END 
LOADS  AND  MOMENTS  WITH  ROTATIONAL  RESTRAINT 


Flat  Bar 

Mode 

m 

R 

c 

Tripping  Stress  [ (0{ 
Coefficient 

)  x  104]/E 

i  ere 

Finite  Element 

Eqn.  (22) 

i 

0 

6.44 

6.42 

2 

0 

7.58 

7.51 

3 


9.33 


0.427 


17.6 

19.3 

10.5 

10.7 

10.9 

10.8 

41.2 

70.6 

16.9 

23.6 

14.1 

16.5 

Tripping  Moment  M  (d- 

-h+  \  t)  x  104  /  [IE] 

Coefficient 

Finite  Element 

Eqn.  (36) 

9.04 

8.99 

10.66 

10.52 

24.3 

27.0 

14.6 

15.0 

55.2 

98.9 

22.8 

33.0 

TABLE  4  -  TRIPPING  STRESS  COEFFICIENTS  FOR  TEES  UNDER  END 
LOADS  AND  MOMENTS  WITH  ROTATIONAL  RESTRAINT 


Tee  Stiffener 

Mode 

m 

R 

c 

4 

Tripping  Stress  [  (o  )  *  10  ]/E 

e  ere 

Coefficient 

Finite  Element 

Eqn.  (21) 

i 

0 

26.4 

26.3 

2 

0 

65.9 

72.6 

1 

0.826 

45.3 

48.0 

2 

0.0748 

73.3 

78.0 

1 

1.65 

56.4 

69.7 

2 

0.150 

77.7 

83.4 

1 

4.96 

75.5 

156.6 

2 

0.449 

85.6 

105.2 

1 

8.26 

82.7 

243.5 

2 

0.748 

88.6 

126.9 

Mode 

m 

R 

Tripping  Moment  M  (d 

Coefficient  Cre 

,-h+  |t)x  104j  j  [ IE] 

Finite  Element 

Eqn.  (35) 

1 

0 

29.4 

29.0 

2 

0 

75.1 

80.1 

1 

0.826 

49.7 

53.0 

2 

0.0748 

82.1 

86.0 

1 

1.65 

61.4 

76.9 

2 

0.150 

86.3 

92.0 

TABLE  5  -  TRIPPING  LOAD  COEFFICIENTS  UNDER  UNIFORM  LATERAL  LOAD 
WITH  ROTATIONAL  RESTRAINT 


Flat  Bar 

2  l  t  l  \ 

a  ^l-li+  ,  l  ) 

Tripping  Load  100  <j 

Mode 

K 

H  .  * 

Cocf  t  lc  lent 

Finite  K lenient 

Kqn.  (47) 

i 

0 

SS 

-1.13 

-  1.12 

1  , 

0 

SS 

-  1.20 

-  1.18 

0 

c 

+  1.76 

+  1 . 46 

Hi 

0 

c 

+  1.  18 

+  3.40 

2 

0.427 

SS 

-  1.49 

-  1.53 

i? 

o.  m 

SS 

-  1.51 

-  1.61 

i 

2.0 

(' 

+  1.15 

+  4.  14 

2 

0.427 

c 

+  1. 69 

+  4.41 

t? 

0.  764 

SS 

-  1 . 8‘l 

-  2.21 

B 

0.  118 

SS 

-  1.41 

-  2.  10 

Eg  | 

2.  14 

c 

+  4.44 

+  8.  J5 

19 

0.  764 

i: 

+  4.  Id 

+  8.18 

Tee  Stiffener 

lvlppliiK  Loud  100  4  ,ri,  Y  l 

Mode  ^ 

Goef 1 lc lent 

Finite  F.lemeut  1 

Kqn.  (46) 

Hi 

0 

SS 

-  4.  59 

-  4.20 

i  ! 

0 

SS 

-  9.0} 

-10.04 

0 

c: 

+  6 .  54 

+  7.41 

1 1 

0 

c 

+17.91 

+40. 7b 

n 

0.826 

SS 

-  (l.  70 

-  8.  75 

9 

0.0748 

SS 

-  4.65 

-10. 74 

I 

0.826 

c 

+  10.  7 

+  1  1.  1 

n 

0.0748 

c 

+20.  1 

+42.  7 

B 

8.2b 

SS 

-10.  1 

-70.7 

i 

0.  748 

SS 

-11.0 

-  1  (i  .  8 

n 

8.26 

c: 

+17.3 

+46.2 

B 

0.  748 

c 

+24.4 

+  48.  1 

*B.C.  •  Boundary  Conditions,  SS  •  Simple  Support, 

(!  -  Clamped. 

7  -  Indicates  uncertainty  as  to  mode  number. 
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significant  source  of  error  for  the  case  of  Lateral  loading,  particularly 
positive  pressure,  is  the  distribution  of  out-of-plane  displacements  in 
the  x-directlon  as  defined  by  Equation  (40).  The  presence  of  this  source 
of  error  is  most  detectable  in  these  examples  for  flat  bars  tuivlng  zero 
rotational  restraint  about  their  line  of  attachment  since,  in  this 
situation,  little  if  any  distortion  of  the  web  occurs  and,  therefore,  the 
Lack  of  clo  ie  agreement  with  the  finite  element  results  cannot  be  ascribed 
to  the  neglect  of  web  deformation. 

Despite  these  discrepancies  noted,  Equation  (40)  does  a  good  Job, 
qualitatively  at  least,  in  describing  the  tripping  phenomena  for  modes 
m  =  l  and  m  »  2  when  rotational  restraint  is  nonexistent  or  very  small. 

As  rotational  restraint  Increases,  the  tripping  patterns  become  in¬ 
creasingly  complex  and  soon  this  complexity  reaches  a  stage  where  it  is 
exceedingly  difficult  to  associate  a  specific  mode  number  with  a  particular 
buckling  pattern,  as  the  question  marks  appearing  in  Table  5  indicate. 

This  is  due  to  the  concentration  of  tripping  deformations  in  the  regions 
of  the  stiffener  undergoing  compressive  loading,  with  strongly  attenuated 
or  no  deformations  at  all  in  those  regions  experiencing  tension.  Thus, 
tripping  deformations  under  negative  lateral  loading  will  tend  to  be 
concentrated  in  the  midsection  of  the  stiffener  while  those  associated 
with  positive  lateral  loading  will  be  concentrated  at  the  two  ends.  In 
fact,  for  positive  lateral  loading  and  practical  geometries  (assuming  there 
is  sufficient  vertical  fixity  to  produce  regions  of  compressive  stresses) 
there  are  effectively  only  two  buckling  modes,  one  symmetrical  and  the 
other  ant isymmetr ical  about  the  stiffener  midpoint.  As  the  degree  of 
rotational  restraint  increases,  the  deformations  associated  with  both 
modes  become  more  and  more  concentrated  towards  the  ends  until  for  large 
degrees  of  restraint  (particularly  for  long  stiffeners)  the  buckling 
patterns  may  be  completely  localized  in  the  compression  regions.  As  the 
buckling  patterns  become  more  localized,  tin'  corresponding  tripping  loads 
approach  a  common  value.  Consequently,  for  long  stiffeners  (large  a/d) 
with  moderate  to  high  rotational  restraint,  the  tripping  load  for  both 
modes  are  effectively  identical  with  the  mode  shapes  essentially  differing 
only  with  regard  to  their  conditions  of  symmetry. 
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The  rapidly  deteriorating  accuracy  of  the  tripping  solutions  described 
in  this  section  as  rotational  restraint  increases,  clearly  indicates,  even 
though  based  on  very  limited  geometries  and  data,  that  these  solutions  are 
at  best  suitable  for  problems  with  zero  or  very  limited  amounts  of  such 
restraint.  Because  of  the  number  of  variables  involved,  parametric  studies 
would  be  required  to  effectively  define  the  degree  of  rotational  restraint 
below  which  an  acceptable  level  of  accuracy  could  be  guaranteed.  Since 
the  use  of  finite  element  programs  in  such  studies  can  present  significant 
economic  stumbling  blocks,  an  alternate  approach  was  adopted  which  led  to 
the  development  of  approximate  analytical  solutions  which  take  into  account 
web  deformations  and  are  suitable  for  either  hand  calculation  or  programming 
on  desk-top  computers.  These  solutions  are  described  and  evaluated  in  the 
next  section.  Unfortunately  these  solutions  are  available  only  for  the 
cases  of  axial  end  load  and  constant  moment,  and  not  for  lateral  loading. 
Consequently,  semiempir ical  solutions  will  probably  have  to  be  developed 
for  tripping  under  lateral  load,  making  use  of  the  analytical  work  pre¬ 
sented  here,  finite  element  results,  and  any  other  relevant  data  available, 
at  least  until  an  improved  analytical  solution  for  this  particular  loading 
case  is  developed. 

EFFECTS  OF  WEB  DEFORMATIONS 

In  the  previous  development  the  critical  assumption  was  made  that  the 
centerline  of  the  web  at  any  transverse  cross  section  would  not  curve, 
l.e.,  for  small  deflections,  v®  zfl.  As  the  selected  finite  element  so¬ 
lutions  have  demonstrated,  however,  this  assumption  in  certain  circum¬ 
stances  can  lead  to  buckling  predictions  which  are  considerably  higher  than 
the  true  values.  Since  the  proportions  of  the  two  stiffener  examples 
studied  are  not  untypical  of  ship-type  structures,  it  is  clear  that  so¬ 
lutions  which  take  into  account  web  deformations  would  be  most  valuable. 

If  the  web  is  allowed  to  curve,  then  for  small  deformations  the  dis¬ 
placement  of  any  point  on  the  centerline  of  the  web  can  bo  described  in 
terms  of  a  translation  v  and  a  rotation  3 .  If  a  deformation  function  is 
assumed  for  the  web,  then  the  displacement  of  any  point  on  the  web 
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centerline  can  be  described  in  terms  of  the  v  and  0  of  a  single  chosen 
point,  in  this  case  the  shear  center,  as  shown  in  Figure  7.  (In  the 
following  development  the  location  of  the  shear  center  is  assumed  to  be 
at  the  intersection  of  the  web  and  flange  centerlines,  i.e.,  s  ®d  ,  and 
although  only  approximately  true,  the  shear  center  is  indeed  very  close  to 
this  point  for  stiffeners  with  flanges  of  typical  proportions  and  the 
consequences  of  this  approximation  are  insignificant.)  A  logical  deforma¬ 
tion  function  is  tliat  of  a  beam,  cantilevered  at  one  end  (the  flange  end) 


Figure  7  -  Characterization  of  Web  Deformations 


33 


and  loaded  by  a  concentrated  force  and  moment  at  the  other  (corresponding 
to  the  restraining  force  and  moment  generated  by  the  plating  to  which  the 
stiffener  is  attached).  In  terms  of  the  v  and  0  displacements  of  the 
shear  center  this  assumption  gives  the  following  expression  for  the  lateral 
displacement  of  any  point  on  the  stiffener  web, 


The  strain  energy  stored  in  an  element  of  the  web  due  to  its  deforma¬ 
tion  as  described  by  Equation  (57)  is  given  by 


(1+3R2)  v2  -  2(1+R2)d  v0 

c 


+ 


2 


dx 


(59) 


The  total  strain  energy  stored  in  the  structure  in  its  buckled  state 
is  thus 
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ei  v2  +  Ere2  +  gj32 

Z  XX  XX  X 


+  Cy2  +  -y-  j(l+3R2)  V2  -  2(1+R2)dcve 
**c 


(l  +  i  R2)  d262  } 


where  y  is  the  rotation  of  the  stiffener  base  (and  thus  the  elastic 
spring)  as  shown  in  Figure  7.  It  is  not  difficult  to  show  that 


(3  S'  -6) 


and,  therefore,  V  completely  in  terms  of  v  and  $  is  given  by 


^  “  d 

v  -  f  ei  v2  +  Ere2  +  cje2  +  -4  r(i-r)  (9v2-6d  ve*d2e2) 

2  J  z  xx  xx  x  ,3  c  c 

0  |_  c 


+  ^y  |(1+3R2)  V2  -  2(1+R2)  dcve  +  (l+  J  R2)  d2e2j 
dc  - 


If  the  stiffener  is  loaded  by  an  axial  inplane  force,  the  potential 
energy  of  the  loading  in  this  case  is  the  negative  of  the  work  done  by 
that  loading,  namely 


uw  -  -  w  -  -  J*J*°X  <S(y,z)dydx 
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0 


(68a) 


<j>  -  4-  $  +  ~ — r  * 

v  dx  vx  d7  Vxx 


*8-sr*s  +7T»e  -° 

x  dx  xx 


(68b) 


where  the  subscript  indicates  a  partial  derivative  with  respect  to  that 
parameter. 

Carrying  out  the  operations  indicated  by  the  above  expressions  and 
assuming  solutions  for  v  and  0  of  the  form 


.  mTTX  ,  a  a  .  mTTX 

v  ■  v  sin  -  and  0  -  0  sin  - 

o  a  o  a 


leads  to  a  set  of  two  homogeneous  equations  for  the  constants  vq  and  0q. 
Nonvanishing  solutions  for  vq  and  0q  can  only  exist  if  the  determinant 
of  the  coefficients  of  these  equations  is  zero.  Computing  this  determinant 
results  in  a  quadratic  equation  for  the  critical  buckling  stress,  (op)  e, 


«Vcr.  <k2kCk6>  +  <*1 VWW  +  <klk3-kI>  '  0  (70) 


where  the  k^'s  are  given  by 


2  D 


kl-EIz(?)  +3^(ti)  (1+3R) 


(71a) 


/ 18  19  3  2\ 

k.  -  -  A  +  d  t  (41  +  R-  -Tj-x  R  ) 
2  s  c  w  \ 35  140  140  / 


(71b) 


k3  ■ GJ  +  Er(r)2  * irHs;)2  (1+iR) 


(71c) 
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k 


5 


k  a  k  — 
*6  K6 


d3t 
c  w 


d  -  h  +  ^t 
c  2 


/!13_ 

l  3360 


140  R+  224 


r2\ 

224  K  / 


(73e) 


(73f) 


Although  the  logic  is  clearly  less  soundly  based,  the  web  deformation 
function  described  by  Equation  (57)  can  also  be  applied  to  the  problem  of 
the  tripping  of  flat  bars.  In  this  case,  v  and  3  represent  not  the  dis¬ 
placement  of  the  shear  center,  but  rather  that  of  the  outer  extremity  of 
the  flat  bar.  The  strain  energy  is  given  by 


it  a 

KJJ 


(l<vv)x*+(vv,)2Zl 


0  0 


J 


R(l-R)  (9v2-6dvB4d2B2)dx 

A  J 


(74) 


while  the  potential  energy  is  identical  to  Equation  (64)  with  the 
parameter  d^  replaced  by  the  flat  bar  depth  d  and  I  replaced  by  the  polar 
moment  of  inertia  of  the  flat  bar  about  its  outer  extremity  j  d3t  )  . 

The  procedure  followed  to  determine  the  critical  stress  is  the  same  as 
that  described  previously  and  results  in  a  quadratic  equation  identical  to 
Equation  (70).  However,  in  this  case,  the  kj  * s  are  defined  as  follows: 


ki  ■  f  (a )  \h  <l«*> + 1 (r)  f*-™'2) 


2  /  md  \  *  /  17 
+  7T  (r)  V  35  ~ 


140  R+  140  R  ) 


(75a) 
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k  *  -  d  t 
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— 2(r)4(x§s »2>; 

W  )tM  -  fe  R+  450  r2+  aio  «"3»«2)  <t„/J>2 


k5--D„(s)  [7  <1+R>+(;r)  (I*'"!1' 


(75b) 


(75c) 


(75d) 


+  l«2)«2(fi)4(3f-flo-llok2)]  <"•> 

k6  ■  +  d\  h  -  m  R+  TTo  r2+  Ho  (2-2R+R2)  <t„/i)2]  (75f) 


(75f) 


For  a  constant  moment.  Equation  (72)  is  applicable.  In  this  case, 
however,  the  k^  functions  are  related  to  the  k^'s  defined  as  Equations 
(75a-75f)  as  follows: 


(=*) 


(76a) 
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+  4 1 L 


39_  _3_  2 
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(76c) 
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_29 _ n.  R+  _J_  r2 

3360  1680  672 


48  ( 10  ~  15  R+  10  R  )  (tw/d) 


(76d) 


(d-h+  |  t 


(76e) 
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k6  ‘ 


d  -  h  +  -z  t 


113 _ 3_  JL_  2 

3360  140  224 


+  48  (  5  "  10  R+  10  R  )  (tw/d) 


(76f) 


While  these  solutions  for  flat  bars  are  approximate,  they  have  the 
advantage  of  not  requiring  an  iterative  solution  (if  be  and  C  are  known). 
An  "exact"  solution  for  end  loading  exists  based  on  thin  plate  theory. 
The  critical  stress  is  determined  from  the  solution  of  the  equation 


[c^  cos  {3d  +  c.,  cosh  ad  +  c^2  sln''  Ic,,  $  cos  0d 


+  c^  0  cosh  ad]  -  [-C2  P  sin  0d  +  c^  a  sinh  ad 


+  c^c^  a  cosh  ad]  |c^  sin  0d  +  cn  ^  sinh  adj  ■  0 


(77a) 


where,  in  this  case 
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2  /  ran  S  ‘ 
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(77b) 


The  obvious  problem  with  this  approach  is  that  an  iterative  solution  of  a 
transcendental  function  is  required,  thus  effectively  eliminating  it  as  a 
practical  tool  for  hand  calculation.  Solution  by  computer  is  quite 
straightforward,  nevertheless  the  necessity  of  an  iterative  solution  may 
still  be  a  serious  drawback  particularly  in  computer  oriented  design 
synthesis  programs.  Several  levels  of  iterations  may  be  required  in  such 
programs,  thus  resulting  in  the  need  for  a  single  critical  stress  to  be 
evaluated  literally  hundreds,  or  even  thousands  of  times.  In  contrast, 
the  solutions  represented  by  Equations  (70)  and  (72)  are  not  inherently 
iterative  by  nature,  and  although  admittedly  cumbersome,  are  suitable  for 
hand  calculation.  Even  when  the  assumption  that  certain  parameters, 
usually  bg  and/or  C,  are  load  dependent  requires  an  iterative  solution, 
the  manual  use  of  these  expressions  is  possible. 
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While  solutions  which  take  into  account  the  influence  of  web 
deformations  are  available  for  the  cases  of  axial  end  load  and  constant 
moment,  unfortunately  such  solutions  do  not  presently  exist  for  the  case 
of  lateral  loading  as  previously  indicated.  The  difficulty  is  primarily 
due  to  the  variation  of  with  x  and  the  inadequacy  of  a  single  sine 
function  to  represent  the  buckled  shape.  The  lack  of  a  solution  for  this 
particular  loading  case  is  particularly  bothersome  since  the  limited  data 
presented  in  the  previous  section  suggests  a  significant  influence  on 
buckling  strength  by  web  deformations.  This  problem  is  presently  being 
investigated. 

Using  the  same  configurations  previously  discussed  (see  Figure  6),  a 
series  of  comparisons  were  made  between  the  solutions  which  include  web 
deformations  and  those  which  do  not.  The  results  of  these  comparisons  are 
presented  and  discussed  in  this  section. 

Tables  6  through  9  present  buckling  coefficients  for  both  tee  and  flat 
bar  stiffeners  of  constant  length  with  varying  degrees  of  rotational 
restrain  about  their  toes.  These  tables  clearly  reemphas  ze  behavior  that 
was  noted  in  the  previous  section,  namely  the  increasing  significance  of 
web  deformations  on  buckling  strength  as  rotational  restraint  increases. 
These  tables  also  indicate  excellent  agreement  between  the  finite  element 
solutions  (the  numbers  on  the  tables  in  parentheses)  and  the  approximate 
solutions  represented  by  Equations  (70)  and  (72)  for  all  the  degrees  of 
restraint  considered.  With  regard  to  the  tee  stiffener  configuration, 
however,  »t  is  quite  noticeable  that  for  small  or  zero  rotational 
restraint.  Equations  (70)  and  (72)  predict  somewhat  lower  (by  several 
percent)  buckling  coefficients  than  do  the  finite  element  solutions.  It 
would  appear  that  this  may  be  primarily  due  to  the  fact  that  in  the 
finite  element  model  the  flange  of  the  tee  is  modeled  bv  a  bar,  which  is  a 
one-dimensional  line  element.  Consequently,  on  computing  the  work  W  of 
the  external  forces,  the  finite  elemeitt  method  is  unable  to  include  the 
effects  on  the  parameters  I  and  S  (see  Equations  (20)  and  (34))  of  the 
finite  width  and  thickness  of  the  flange.  As  the  degree  of  rotational  re¬ 
straint  increases,  the  data  strongly  suggest  that  the  Influence  of  this 
approximation  decreases  in  significance. 
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TABLE  6  -  TRIPPINi:  STRESS  COEFFICIENTS  FOR  TEES  UNDER  END  LOADS 
WITH  ROTATIONAL  RESTRAINT 


Toe  Silt  toner 


Mode 

in 

R 

c 

Tripping  Stress 
Coef f ic lent 

|(o  )  x  104]/E 

e  ere 

Eqn. 

(70)* 

Eqn.  (21) 

i 

0 

24.3 

(26.4) 

26.3 

2 

0 

65.5 

(65.9) 

72.6 

l 

0.826 

43.5 

(45.3) 

48.0 

•> 

0.075 

7  J.l 

(73.3) 

78.0 

i 

1.65 

54.5 

(56.4) 

69.  7 

) 

0.150 

77.6 

(77.7) 

83.4 

i 

4.  9h 

73.4 

(75.5) 

156.6 

•) 

0.449 

85.4 

(85.6) 

105.2 

i 

8.26 

80.4 

(82.7) 

243.5 

) 

0.  748 

88.3 

(88.6) 

126.9 

i 

15.0 

86 . 3 

420.7 

2 

1 . 36 

90.9 

171.2 

♦Figures  in  parent  hoses  arc  finite  element 
values. 


Flat  Bar 

Mode 

m 

R 

c 

Tripping  Stress  [ (o 
Coefficient 

)  x  io4]/e 

ere 

Eqn.  (70)* 

Eqn.  (22) 

i 

0 

6.42  (  6.44) 

6.42 

2 

0 

7.51  (  7.58) 

7.51 

1 

1.00 

12.44 

12.84 

2 

0.214 

9.10 

9.12 

1 

2.00 

17.6  (17.6) 

19.3 

2 

0.427 

10.5  (10.5) 

10.7 

1 

5.00 

29.3 

38.5 

2 

1.07 

13.6 

15.5 

1 

10.0 

41.5  (41.2) 

70.6 

2 

2.14 

17.0  (16.9) 

23.6 

1 

15.0 

49.1 

102.7 

2 

3.70 

_ 

19.1 

31.6 

*Flgures  in  parentheses  are  finite 
values. 

element 

TABLE  9  -  TRIPPING  LOAD  COEFFICIENTS  FOR  FLAT  BARS  UNDER  END  MOMENTS 

WITH  ROTATIONAL  RESTRAINT 


Flat  Bar 


Mode 

m 

R 

c 

Tripping  Moment 
Coef  f ic lent 

(m  (d 

L  ere  \ 

-h+  |  t)  x  io4j  1  [IE] 

Eqn. 

(72)* 

Eqn.  (36) 

i 

0 

8.99 

(  9.04) 

8.99 

0 

Am 

0 

10.50 

(10.66) 

10.52 

1 

1.00 

17.2 

18.0 

2 

0.214 

12.6 

12.8 

i 

2.00 

24.2 

(24.3) 

27.0 

2 

0.427 

14.4 

(14.6) 

15.0 

i 

5.00 

39.6 

54.0 

2 

1.07 

18.5 

21.8 

1 

10.0 

55.3 

(55.2) 

98.9 

2 

2.14 

22.7 

(22.8) 

33.0 

1 

15.0 

64.8 

143.9 

1 

3.20 

25.3 

44.2 

*Figures  in  parentheses  are  finite  element  values. 


It  was  noted  In  the  previous  section  that  not  only  degree  of  ro¬ 
tational  restraint  but  also  stiffener  proportions  can  influence  the 
amount  of  web  deformation  which  will  be  present.  An  example  of  this  is 
illustrated  in  Table  10  where  tripping  stresses  have  been  calculated  for 
the  tee  stiffener  of  Figure  5  for  varying  lengths.  It  may  be  noted  that 
as  the  length-to-depth  ratio  decreases  there  is  an  increasingly  large 
discrepancy  between  the  values  of  the  tripping  stresses  predicted  by 
Equations  (21)  and  (70),  indicating  an  increasing  role  of  web  deformation. 
As  has  been  noted,  tripping  for  very  short  lengths  is  characterized  by 
deformations  consisting  mainly  of  buckling  of  the  web  and  twisting  of  the 
flange  about  its  centerline.  This  is  clearly  a  local  buckling  phenomena 
and  its  presence  is  indicated  in  Table  10  by  these  cases  for  which  large 
discrepancies  exist  between  the  finite  element  results  and  those  of  Equa¬ 
tion  (70).  These  discrepancies  occur  because  Equation  (57),  while  an 
excellent  approximation  to  the  deformations  which  occur  in  the  web  during 
tripping,  is  a  somewhat  less  than  ideal  approximation  for  local  buckling  of 
the  web.  Thus,  while  Equation  (70)  has  been  demonstrated  to  be  reliable 
for  predicting  the  primary  tripping  of  tee  stiffeners,  it  seems  clear  from 
this  example  that  the  possibility  of  local  web  buckling  must  be  investi¬ 
gated  by  other  means,  such  as  plate  theory. 

For  tee  stiffeners,  other  geometrical  parameters  besides  a/dc  can 
influence  the  amount  of  web  deformation  present.  In  a  qualitative  sense, 
the  most  significant  factor  involved  seems  to  be  the  relative  stiffnesses 
of  flange  versus  web,  since  a  stiff  flange  will  tend  to  resist  deformation, 
therefore,  forcing  more  into  the  web.  Many  details  of  geometry  can  in¬ 
fluence  this  ratio  and  thus,  no  further  attempts  will  be  made  to  examine 
them  here.  However,  Table  11  illustrates  the  influence  of  length  on 
degree  of  web  deformation  when  there  is  no  flange.  The  absence  of  any 
significant  influence  of  length  on  tripping  in  this  case  in  essence 
supports  the  hypotheses  of  the  importance  of  flange  versus  web  stiffness 
relative  to  web  deformation. 
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TABLE  10  -  TRIPPING  STRESS  COEFFICIENTS  FOR  TEES  UNDER  END  LOADS  WITH 
VARYING  LENGTH/DEPTH  RATIOS  (NO  ROTATIONAL  RESTRAINT) 


Tee  Stiffener 

4 

Tripping  Stress  [(0p)rrt,  x  10  1/E 

Mode 

a 

Coefficient 

m 

c 

Eqn.  (70)* 

Eqn.  (21) 

1.99 

320.4  (193.1) 

449.8 

2.66 

217.4  (188.2) 

257.7 

3.99 

107.8  (105.2) 

120.6 

4.06 

104.7  (102.5) 

117.0 

5.32 

65.5  (  65.8) 

72.6 

7.48 

38.3 

42.1 

8.11 

34.2  (  35.8) 

37.4 

10.64 

24.3  (  26.4) 

26.3 

14.96 

17.6 

18.7 

16.62 

16.3 

17.2 

— 

33.24 

12.2 

12.4 

*Figures  in 

parentheses  are  finite  element 

values. 
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TABLE  11  -  TRIPPING  STRESS  COEFFICIENTS  FOR  FLAT  BARS  UNDER  END 
LOADS  WITH  VARYING  LENGTH/DEPTH  RATIOS  (NO 
ROTATIONAL  RESTRAINT) 


Flat  Bar 

4 

Tripping  Stress  [  (o  )  x  io  )/E 

Mode 

a 

Coefficient 

m 

md 

Eqn.  (70)* 

Eqn.  (22) 

1 

20.3  (20.0) 

20.3 

2 

9.61  (  9.51) 

9.61 

3.12 

7.51  (  7.47) 

7.51 

4 

6.95  (  6.92) 

6.95 

5 

6.62  (  6.61) 

6.62 

6.25 

6.42  (  6.41) 

6.42 

m 

6.31  (  6.30) 

6.31 

m 

6.20  (  6.20) 

6.20 

15.0 

6.12  (  6.12) 

6.12 

30.0 

6.07  (  6.07) 

6.07 

— 

♦Values  in  parentheses  from  Equation  (77). 
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TRIPPING  UNDER  COMBINED  LOADING 

The  treatment  of  tripping  to  this  point  has  considered  both  inplane 
and  lateral  loadings,  but  only  applied  individually.  In  many  practical 
situations  such  loads  will  occur  simultaneously  which  naturally  leads  to 
the  subject  of  tripping  under  combined  loading. 

One  approach  to  this  problem  is  to  follow  the  same  procedure  used  for 
the  individual  loading  cases,  that  is,  to  equate  the  strain  energy 
associated  with  the  assumed  deformation  function  to  the  work  done  by  the 
external  loads.  In  this  case,  however,  the  work  will  include  contributions 
from  all  the  simultaneously  occurring  loads.  One  stress  (or  load)  is  then 
solved  for  (usually  the  predominant  one)  and  thus  becomes  the  tripping 
stress  (or  load)  "in  the  presence  of"  the  other  applied  loads. 

As  an  example,  consider  the  case  of  a  tee  stiffener  loaded  by  both  an 
axial  end  load  and  a  uniform  lateral  load.  Since  a  lateral  load  is  present, 
the  displacement  function  represented  by  Equation  (AO)  must  be  used.  The 
strain  energy  associated  with  this  displacement  function  has  previously 
been  computed,  but  is  repeated  here, 


V  =  7  a  B2  (-)  GJ  H  (K)  (78) 

A  o  \a  /  m 

for  convenience.  The  work  done  by  the  lateral  load  has  also  been  previously 
determined, 


W 

q 


-  Is  ( i )  so  "2  *  s  F«(l° 


(79) 


but  here  it  is  represented  by  the  symbol  to  indicate  it  is  the  component 
due  to  the  lateral  load.  Equation  (20)  gives  the  work  done  by  the  axial 
end  load  for  the  assumed  displacement  function  given  by  Equation  (12); 
for  displacement  Equation  (AO)  the  appropriate  work  expression  becomes 


“P  -K  »lil)  »  yA2-K-«>2  <k-»2] 


(80) 
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The  total  work  Is  the  sum  of  the  two  components. 


W  -  W_  +  W 

P  q 

■  i  °e  (;)  a  Ipl»V+(mt2)2  (K-l)2) 

~h  (!)bo  1,2  a  S  VK>  <81> 


Equating  the  total  work  and  the  strain  energy  defines  the  condition 
of  loading  at  which  tripping  will  occur.  If  the  axial  end  load  is  assumed 
to  predominate,  then  carrying  out  the  above  procedure  and  solving  for  0 
produces  the  following: 


GJ  Hm<K)  +  12  \l)  S  Fm(K) 
Iplm2K2+(m+2)2  (K-l)2] 


(82) 


As  in  the  case  of  lateral  load  alone,  the  parameter  K  must  be  selected  so 

that  the  resulting  value  of  (o  )  is  the  minimum.  Since  both  the 

e  ere 

numerator  and  denominator  of  Equation  (82)  can  be  represented  as  quadratic 
functions  of  K,  the  procedure  previously  described  for  computing  the 
appropriate  value  of  K  can  also  be  applied  here. 

A  second  approach  for  computing  tripping  under  combined  loading  in¬ 
volves  the  application  of  the  so-called  "interaction"  formula.  In  this 
approach,  for  the  type  of  loads  considered  here,  tripping  is  assumed  to 
occur  when  the  following  condition  is  satisfied: 
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are  the  failure  "loads" 


In  this  expression  (o  )  ,  M  ,  and  q 

e  ere  ere  xre 

corresponding  to  separate  application  of  the  three  types  of  loads* 
considered  (and  all  for  the  same  mode  number  even  though  minimum  failure 
loads  for  the  different  loads  may  occur  for  different  modes).  The 
indices  a,  3,  and  Y»  while  often  suggested  by  theoretical  considerations, 
are  usually  selected  empirically  on  the  basis  of  experimental  and/or 
numerical  results. 

Limited  calculations  with  expressions  such  as  Equation  (82),  an 
example  of  which  is  illustrated  in  Figure  8  for  the  "tee"  of  Figure  6, 
suggest  that  for  design  purposes  a  linear  interaction,  as  follows. 


0  0.2  0.4  0.0  0.0  1.0 

q 


qcr« 

Figure  8  -  Interaction  Diagram  for  Combined  Axial  and 

Lateral  Loads 


*The  distinction  between  M  and  the  moments  associated  with  q  is  the 
following.  The  moments  associated  with  q  arise  due  to  local  bending  of  the 
plate-beam  between  supporting  transverses,  whereas  M  is  the  moment  compo¬ 
nent,  approximately  constant  between  adjacent  transverses,  due  to  overall 
grillage  bending. 
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M 


1 


(84) 


<°e>cre 


ere 


+  -S-- 

^cre 


Is  the  appropriate  one  to  use.  Although  a  bit  conservative,  this  selection 
seems  warranted  until  further  evidence,  either  analytical  or  experimental, 
indicates  a  more  realistic  choice. 

Generally  speaking,  the  worst  loading  condition  which  can  occur  is 
one  in  which  all  the  loads  present  are  directed  so  as  to  have  a  de¬ 
stabilizing  effect  on  the  structure,  i.e.,  the  signs  of  all  the  ratios  in 
Equation  (84)  are  positive.  Certain  loads,  however,  may  effectively 
contribute  to  structural  instability  only  when  they  act  in  a  preferred 
direction,  such  as  a  compressive  stress  0£  or  a  positive  moment  M,  for 
example.  When  such  loads  act  in  the  opposite  direction  they  will  produce 
a  stabilizing  effect  in  the  structure  which  will  be  indicated  by  a  negative 
sign  for  the  appropriate  ratio  in  Equation  (84).  In  such  cases,  the  worst 
loading  condition  will  likely  be  one  in  which  only  some  of  the  loading 
components  are  present.  The  choice  of  what  loading  components  to  include 
in  design  calculations  and  whether  to  take  into  account  stabilizing 
effects  is  a  decision  the  designer  must  make  on  a  case-by-case  basis 
relative  to  his  knowledge  of  what  loading  combinations  can  and  are  likely 
to  occur. 

Relative  to  the  choice  of  which  approach  to  adopt  in  treating  tripping 
under  combined  loading,  the  interaction  formula  approach  offers  certain 
practical  advantages.  Empirical  or  semiempirical  solutions,  as  well  as 
those  of  a  more  vigorous  analytical  nature,  are  readily  adaptable  to  an 
interaction  formulation.  This  is  particularly  significant  in  this  case 
since  tripping  expressions,  which  take  web  deformations  into  account,  are 
not  presently  available  for  lateral  loading.  Semiempirical  solutions, 
therefore,  will  have  to  be  employed,  at  least  on  a  temporary  basis.  The 
form  of  the  interaction  expression,  through  the  indices  ex,  3,  and  y>  also 
can  readily  accommodate  the  results  of  experimental  investigations,  if  and 
when  such  data  are  available  (and  regardless  of  what  values  any  analytical 
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theories  say  such  indices  should  assume).  This  great  flexibility  and  the 
ease  of  application  that  the  interaction  formulation  offers  make  this 
approach  a  very  reasonable  and  practical  choice  for  design  applications. 

ESTIMATION  OF  PLATE  ROTATIONAL  RESTRAINT 
The  numerous  examples  presented  more  than  adequately  illustrate  the 
significant  increases  in  tripping  strength  that  rotational  restraint  along 
the  toe  of  the  stiffener  provides.  To  take  advantage  of  this  added 
strength,  it  is  necessary  to  be  able  to  estimate  the  degree  of  rotational 
restraint  that  the  plating,  to  which  the  stiffener  is  attached,  provides. 

One  approach  to  this  problem  is  to  apply  Rayleigh’s  principle  as 
previously  described  but  to  include  one  frame  spacing  of  plating  in 
addition  to  the  stiffener.  When  tripping  occurs  under  axial  end  load  the 
plating  is  assumed  to  deform  laterally  according  to  the  expression 

w  =  w  sin  sin  ^  (85) 

o  a  b 

Since  the  junction  of  the  stiffener  web  and  the  plating  remain  at  right 
angles,  this  expression  can  readily  be  written  in  terms  of  the  rotation  B 
of  the  stiffener,  giving 


Computing  the  strain  energy  of  the  plate-stiffener  combination  and  the 
work  done  by  the  external  forces  and  equating  the  two  allows  the  critical 
tripping  stress  to  be  determined  as  follows, 
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GJ  +  E(Iss2fD  (51)  + 


(a  ) 
e  ere 


I  +^-^- 
P  ~  2  ,  xm 
2tt  (o  )  , 
e  pbe 


where  o™bg  is  the  classical  elastic  plate  buckling  stress  for  mode 
number  m. 


■?*•$(&) 2 


and  ^°e)p{,e  *-s  c^e  Pea^  value  of  stress  (at  the  stiffener)  at  elastic 
plate  buckling  corresponding  to  a  nonuniform  distribution  of  stress  in  the 
y-direction.  (When  the  distribution  of  stress  in  the  y-direction  is 
uniform,  (oe>pbe  i®  to  cP^.) 

If  it  is  assumed  that  the  rotational  spring  constant  of  the  plate 
interacts  linearly  with  the  axial  inplane  stress. 


o 

1 - - - 

(o  )m. 

e  pbe 


where  Cq  is  the  spring  constant  of  the  unloaded  plate,  then  this  relation¬ 
ship  can  be  substituted  into  Equation  (21)  and  a  modified  expression  for 

(o  )  determined.  This  gives 
'  e  ere  ° 
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An  examination  of  Equations  (87)  and  (90)  indicates  that  they  are 
identical  provided  that 


C 
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2 


(91) 


Solving  the  above  expression  for  Cq  and  making  use  of  Equation  (88)  pro¬ 
vides  the  following 


C 
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1  t/d 

2  b 


2 


(92) 


which,  in  combination  with  Equation  (89),  is  the  relationship  being  sought. 

For  practical  applications,  certain  modifications  of  the  above  ex¬ 
pressions  are  desirable  at  least  until  numerical  and/or  experimental 
validation  of  the  theory  can  be  provided.  One  such  modification  involves 
the  use  of  the  buckling  stress  (cre)pj3e  Equation  (89).  This  stress  will 
not,  in  general,  correspond  to  the  minimum  plate  buckling  stress  since 
the  mode  number  m  used  in  determining  (°e)pj3e  is  the  tripping  mode 
number.  This  mode  number  will  generally  be  lower  than  that  corresponding 
to  the  minimum  plate  buckling  stress,  here  denoted  by  (0e)pbe-  Although 
the  problem  is  being  investigated,  there  remains  no  clear  understanding  of 
the  state  of  restraint  (or  destabilization)  provided  by  the  plating  when 

the  stress  condition  (a  )  ,  <0  <  (o  )m,  holds.  (Also,  if  the  true 

e  pbe  e  e  pbe 

tripping  mode  shape  contains  any  Fourier  components  of  the  mode  correspond¬ 
ing  to  the  minimum  stress  plate  buckling  mode,  it  appears  that  linear 
theory  would  predict  that  the  plate  restraint  would  drop  suddenly  to  zero 
as  the  minimum  plate  buckling  stress  is  approached,  although  probably  not 
in  a  linear  manner.)  Consequently,  it  is  suggested  that  Equation  (89)  be 
modified  to  the  following: 
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(91) 


Since  a  conservative  estimate  of  effective  wlilth  will  lead  to  an 
optimistically  high  (unconservative)  estimate  of  ((,t,)pj)e  one  further 
modification  Is  called  for,  at  least  tin  til  more  realistic  values  of 
effective  width  are  routinely  applied.  This  involves  the  use  of  the 
classical  elastic  plate  buckling  stress  o  defined  by 


pbe 


4?i*' 1) 
tb2 


(94) 


In  place  of  the  value  (op)  This  leads  to  the  formulation  for  C: 


C  (l- 

°  V  °phe  ' 


C  - 


0  <  0  . 
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a  '  o  . 

e  —  pbe 


II  the  conservatism  of  Equation  ('>*>)  Is  viewed  to  be  excessive,  then  a 

less  conservative  option  is  possible  through  the  use  of  the  average  plate 

stress  0  (not  to  bo  confused  with  the  mean  stress  for  the  cross  section 
a  — 

o  )  In  place  of  o  In  the  above.  However,  it  is  significant  to  note  that 
m  r  e 

if  the  common  practice  of  assuming  a  uniform  distribution  of  primary  stress 

for  stress  levels  below  plate  buckling  is  adopted,  then  -  o  -  o^, 

(o  )  »  0  .  ,  and  Equations  (93)  and  (95)  fire  Identical. 

©pbe  pbe  1 

Because  of  insufficient  knowledge  concerning  the  plat  ing's  behavior, 
one  additional  modification  must  be  considered  witli  regard  to  estimating 


the  appropriate  value  of  C.  The  expression  for  ll  ,  Equation  (92),  was 
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derived  on  the  basis  of  axial  end  loading  and  assumes  that  the  tripping 

mode  shape  was  a  constant  amplitude  sinusoid.  For  axial  end  load  and 

constant  moment  this  is  appropriate,  but  for  lateral  loading  the  more 

complicated  tripping  mode  shapes  which  have  been  noted  and  described  make 

the  application  of  Equation  (92)  for  m  >  i  open  to  question.  Discretion  is 

obviously  called  for  in  this  situation.  Since  C  is  a  monotonically  in- 

o 

creasing  function  of  m,  one  approach  which  can  be  adopted  without  fear  of 
overly  optimistic  predictions  is  to  use  the  minimum  value  of  C^,  correspond¬ 
ing  to  m  ■  1,  namely 


_  1  n‘l) 

Co  "  2  T~ 

for  all  mode  numbers.  While  this  "minimum  restraint"  approach  will 
obviously  result  in  conservative  predictions  in  some  instances,  its 
application  will  still  offer  the  potential  for  significant  increases  in 
tripping  strength  since  the  greatest  relative  increases  are  associated 
with  the  initial  increments  of  rotational  stiffness.  Until  additional 
research  clarifies  these  issues,  this  approach  is  recommended. 

INELASTIC  EFFECTS 

Structures,  in  general,  will  begin  to  exhibit  nonlinear  behavior  at 
stress  levels  below  the  yield  point  due  to  the  presence  of  stress  concen¬ 
trations,  fabrication  distortions,  residual  stresses,  etc.  As  a  result  of 
the  tremendous  complexities  involved  in  treating  this  problem  in  a 
rigorous  manner,  there  is,  at  present,  no  completely  satisfactory  method 
for  considering  the  inelastic  tripping  of  stiffeners  welded  to  continuous 
plating,  particularly  for  application  to  structural  design  calculations. 

It  is  possible,  however,  to  apply  some  simple  "corrections"  to  the  elastic 
solutions  previously  derived  which  are  intended  in  a  gross  manner  to 
approximate  the  effect  of  such  nonlinear  behavior  on  the  tripping  strength 
of  plate-stiffener  combinations. 
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(!) 


S') 


j 

The  approach,  which  is  described  in  more  detail  elsewhere,  involves 
the  application  of  a  parameter  referred  to  as  the  structural  proportional 
limit  ratio  p^.  This  ratio  is  defined  as 

p  -  o  /ov  (97 

r  ps'  Y 

where  0ps  is  the  structural  proportional  limit  stress  and  is  the 

material  tensile  yield  stress.  A  value  of  p  *  0.5  is  typical  for  welded 
3  r 

ships. 

Assuming  that  the  inelastic  tripping  stress  (o  )  is  given  by 


^°e^cr  E  ^°e^cre 


(98) 


and  approximating  the  ratio  E^/E  using  the  Ostenfeld-Bleich  quadratic 
parabolae 


o(oY-o) 

a  (ov-o  ) 
ps  Y  ps 


(99) 


leads  to  the  following  tripping  strength  expressions: 
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(100) 


Although  the  nonuniform  stress  distribution  renders  the  situation 
considerably  more  complex,  this  general  approach  can  also  be  applied  to 
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the  case  of  constant  moment  or  uniform  lateral  pressure.  In  one  such 
treatment,  Equation  (100)  is  still  used,  although,  in  these  cases, 

(ae)cre  is  replaced  by  the  maximum  compressive  bending  stress  corresponding 
to  either  Mcre  or  <lcre  obtained  by  dividing  the  bending  moment  by  the 
section  modulus  Z  of  the  plate  stiffener  section.  The  difficulty  with  this 
approach  is  that  it  inherently  implies  that  some  sort  of  major  structural 
"failure"  will  occur  when  the  maximum  compressive  stress  in  the  outer 
fiber  of  the  stiffener  reaches  the  material  yield  stress.  This  is  clearly 
a  conservative  assumption.  An  alternative  approach  might  be  to  assume 
that  the  inelastic  tripping  moment  Mcr  is  given  by 
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where  (EI)t  might  be  called  a  "tangent  stiffness."  Approximating  the 
ratio  (EI)t/EI  in  a  manner  similar  to  that  of  Equation  (99)  for  Et/E  will 
lead  to  the  following  expressions  for  inelastic  tripping  moment 
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In  these  expressions  is  the  fully  plastic  moment  for  the  plate-stiffener 
combination  and  p^,  while  still  a  structural  proportional  limit  ratio,  now 
represents  the  ratio  of  the  "proportional  limit  moment"  (the  value  of  the 
moment  beyond  which  the  moment-curvature  relationship  of  the  plate- 
stiffener  begins  to  behave  in  a  nonlinear  fashion)  to  the  fully  plastic 
moment.  Although  there  is  little  specific  data  to  indicate  a  value  for  pf 
in  this  case,  a  value  of  p  *0.5  as  recommended  for  use  in  Equation  (100) 
would  seem  to  be  a  reasonable  choice  with  which  to  start. 
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The  degree  to  which  M£r  obtained  from  Equation  (102)  will  differ  from 
that  obtained  from  Equation  (100)  is  primarily  dependent  on  the  shape 
factor  of  the  plate-stiffener  combination.  The  shape  factor  is  the  ratio 
of  the  fully  plastic  moment  of  the  cross  section  to  the  moment  required  to 
cause  initial  yielding.  For  stiffeners  which  have  much  of  their  cross 
sectional  area  concentrated  in  the  flange,  the  shape  factor  will  take  on 
values  not  much  greater  than  1  and,  consequently,  the  differences  resulting 
from  using  the  two  approaches  will  be  minimal.  For  flat  bars,  however, 
the  shape  factor  will  take  on  greater  values  and  the  corresponding  differ¬ 
ences  between  computed  values  of  the  inelastic  tripping  moment  will  be  more 
substantial.  For  example,  the  shape  factors  for  the  cross  sections  illu¬ 
strated  in  Figure  5  are  1.72  for  the  flat  bar  and  1.30  for  the  tee. 

A  similar  approach  can  also  be  applied  in  the  case  of  combined  load¬ 
ing,  with  the  inelastic  "correction"  being  applied  to  the  stresses 
corresponding  to  the  predominant  loading.  Thus,  where  uniform  compressive 
loading  is  predominant.  Equation  (84)  would  be  solved  for  0  (which  in  this 
case  would  represent  the  critical  elastic  tripping  stress  in  the  presence 
of  the  other  loadings)  and  this  stress  then  substituted  in  Equation  (100) 

in  place  of  (0). 

e  ere 

Since  tripping  stresses  must  often  be  calculated  using  an  iterative 
scheme  (when  plate  rotational  restraint  is  included)  there  are  two  possible 
approaches  for  including  the  inelastic  correction.  One  approach  involves 
carrying  out  the  iterations  "completely  elastically,"  including  the  use  of 
Opbe  in  Equation  (95)  for  C,  and  then  applying  the  inelastic  modification 
to  the  converged  solution.  The  alternate  approach  involves  applying  the 
inelastic  modification  within  each  iterative  cycle.  In  this  case  it  would 
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CONCLUSIONS 


This  report  describes  the  application  of  classical  energy  theory  to 
the  problem  of  developing  design  oriented  techniques  for  predicting  the 
lateral  torsional  (tripping)  instability  of  plate-stiffener  combinations. 
Closed  form  solutions  have  been  developed  which  show  excellent  agreement 
with  numerical  results  generated  using  the  finite  element  method  for  axial 
end  loading  and  constant  moment,  but  which  are  less  universally  reliable 
in  treating  plate-stiffeners  under  uniform  lateral  loading. 

The  comparative  solutions  made  with  the  finite  element  analyses  have 
very  clearly  illustrated  the  importance  of  including  the  effects  of  web 
deformations,  particularly  for  the  stiffeners  possessing  substantial 
flanges.  Solutions  which  ignore  these  effects  have  been  shown  to  retain 
acceptable  accuracy  only  when  the  rotational  restraint  present  along  the 
stiffeners  line  of  attachment  to  the  plating  is  nonexistent  or  very 
small.  In  the  course  of  the  development  this  discovery  provided  the  in¬ 
spiration  which  eventually  resulted  in  the  derivation  of  solutions  for 
tripping  under  axial  end  load  and  constant  moment,  including  the  effects 
of  web  deformation.  These  solutions  have  provided  tripping  load  pre¬ 
dictions  which  are  in  excellent  agreement  with  the  corresponding  finite 
element  analyses  for  the  complete  range  of  rotational  restraints  con¬ 
sidered.  They  are  clearly  more  cumbersome  to  use  in  hand  calculations 
than  the  more  simplified  (and  less  accurate)  expressions,  but  their 
complexity  does  not  preclude  their  effective  use  in  manual  calculations, 
particularly  in  view  of  the  increasing  capability  and  sophistication  of 
hand-held  calculators  and  desk  top  mini-computers.  These  solutions  are 
obviously  ideally  suited  for  programming  on  digital  computers  and  thus 
should  be  extremely  effective  when  incorporated  into  large  structural 
synthesis  programs  whose  primary  applications  are  in  the  preliminary 
stages  of  structural  design. 

The  major  problem  area  encountered  in  this  study  is  that  of  estimating 
the  critical  values  for  tripping  under  lateral  loading.  For  the  solutions 
presented,  errors  which  range  approximately  from  0  to  20  percent  can 
reasonably  be  expected  for  the  primary  tripping  mode  (m=l)  with  zero  or 
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extremely  limited  rotational  restraint.  For  larger  degrees  of  rotational 
restraint  and  the  higher  modes  which  are  a  natural  consequence,  the 
corresponding  errors  are  appreciably  larger.  The  inability  at  present  to 
take  web  deformations  into  account  and  to  adequately  model  the  generally 
more  complex  modal  shapes  (apart  from  the  web  deformation  aspect)  are  the 
primary  culprits.  These  difficulties,  in  turn,  are  related  to  the 
functional  dependence  of  the  stress  on  the  x-coordinate.  Until  more 
satisfactory  analytical  solutions  are  developed,  either  semiempirical 
solutions  based  on  the  theory  presented  here  and  numerical  analysis,  or 
finite  element  programs  directed  specifically  at  the  tripping  problem, 
will  have  to  be  utilized. 

Several  other  areas  for  which  further  study  is  required  have  been 
indicated,  or  at  least  alluded  to,  in  the  main  text.  With  regard  to  the 
estimation  of  plate  rotational  restraint,  several  problem  areas  can  be 
identified:  (1)  the  estimation  of  restraint  for  the  more  complex  tripping 

mode  shapes  which  occur  when  lateral  loading  is  present,  (2)  an  improved 
definition  of  the  influence  of  plate  buckling  on  plating  restraint,  and 
(3)  the  effect  of  post  buckling  behavior  of  the  plating  on  stiffener 
tripping.  Inelastic  tripping  behavior  is  another  very  fertile  area  for 
further  research.  It  is  difficult  to  be  specific  in  this  area  since  the 
present  design  oriented  techniques  as  outlined  here  are  extremely  crude 
and  a  careful  examination  of  the  whole  area  is  easily  justified.  Tripping 
of  nonuniform  and  unsymmetrical  stiffeners  are  also  subjects  of  major  and 
practical  interest.  This  "wish  list"  is  not  all-inclusive,  nor  is  it  in¬ 
tended  to  be,  but  it  should  clearly  indicate  that  there  are  certainly 
many  questions  remaining  to  be  answered  with  regard  to  the  tripping 
phenomenon. 

The  ultimate  goal  of  this  study  being  design  applications,  this 
report  would  not  be  complete  without  consideration  of  the  current  Navy 

g 

design  practice  relative  to  stiffener  tripping.  In  this  design  practice 
the  maximum  ratio  of  stiffener  span  to  flange  width  (a/fw),  for  which 
support  is  required  only  at  the  ends,  is  given  by 
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If  che  following  approximations  are  made, 
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then  Equation  (103)  can  easily  be  shown  to  be  identical  to  Equation  (21) 
with  C  =  0,  (a  )  *  a„,  v  =  0.3,  and  m  =  Jl.  The  use  of  this  criteria 

is  designed  to  insure  that  stiffeners  are  proportioned  such  that  tripping 
will  not  occur  at  any  stress  less  than  yield.  This  criteria  is  somewhat 
conservative  since  it  is  based  on  a  theory  for  purely  axial  compression 
(whereas  most  stiffeners  experience  loading  which  is  a  combination  of 
bending  and  axial  compression)  and  it  ignores  the  rotational  restraint  of 
the  plating  (00) .  Factors  which  offset  this  inherent  conservatism  to  a 
small  degree  involve  the  selection  of  m  -  /2  and  the  neglect  of  web 
deformations.  The  value  of  m  chosen  is  intended  to  reflect  the  degree  of 
rotational  restraint  in  the  plane  of  the  web  anticipated  at  the  ends  of 
the  stiffener.  The  neglect  of  web  deformations  in  this  case  is  not 
significant  since,  when  plating  restraint  is  also  neglected,  the  effect  of 
web  deformations  is  in  general  quite  small.  It  should  also  be  noted  that 
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the  criterion  represented  by  Equation  (103)  is  applicable  only  to  tee 
stiffeners  and  for  a  number  of  reasons  cannot  be  appropriately  applied  to 
flat  bar  stiffeners. 

As  a  final  comment,  it  has  become  quite  common  in  reports  of  this 
type  to  decry  the  lack  of  experimental  data  which  are  available  for 
validation  of  the  theories  presented.  This  report  is  no  different.  The 
value  of  good  experimental  data  cannot  be  overestimated  because  any 
predictive  techniques,  no  matter  how  sound  their  theoretical  base,  are 
only  of  value  insofar  as  they  are  capable  of  describing  what  occurs  in  the 
real  world.  The  development  of  confidence  in  the  application  of  new 
technology  in  the  design  process  will  continue  to  be  a  painfully  slow 
process  without  the  validation  only  correlation  with  experimental  data  can 
provide. 
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APPENDIX  A 
SAMPLE  PROBLEMS 

Although  the  application  of  many  of  the  expressions  in  the  main  text 
is  perfectly  straightforward,  the  use  of  sample  problems  can  often  be 
very  enlightening,  particularly  when  the  solution  procedure  requires  a 
number  of  steps. 

Consider  the  plate-tee  stiffener  combination  whose  geometry  is  defined 
in  Figure  6.  A  number  of  geometrical  parameters  require  evaluation  in 
order  for  the  tripping  loads  to  be  determined.  These  parameters  are 
defined  in  Figure  2  and  for  this  example  have  the  following  values: 

I  =  58.79  cm4 
z 

—  =  14.98  cm 
s 

r  =  40.57  cm6 

J  =  9.34  cm4 

I  =  3306  cm4 

P 

ZERO  ROTATIONAL  RESTRAINT 

If  no  rotational  restraint  is  assumed,  i.e.,  C  =  0,  then  a  reasonable 

value  for  the  tripping  stress  for  mode  m  =  1  under  axial  end  loading  may 

6  2 

be  found  using  Equation  (21).  Assuming  values  of  E  =  6.894  *  10  N/cm*’ 
and  v  *  0.3  (which  result  in  a  value  G  =  2.652  x  10  N/cm*"),  this 

expression  gives 

'  ^06  I'2'652  *  10‘')  <9'34) 

+  (6.894  x  106)  ((58. 79)  (14. 98)2  +  40.57)  (tt/100)2] 


The  procedure  for  determining  the  tripping  load  under  lateral  loading 

is  somewhat  more  involved.  First  the  coefficients  of  the  H  (K)  and  F  (K) 

m  in 

functions  must  be  determined.  The  most  convenient  form  of  these  ex¬ 
pressions  is  provided  in  Appendix  B.  Using  these  expressions  the  h^  co¬ 
efficients  can  be  readily  calculated,  giving  for  mode  m  =  1  and  C  =  0  the 
following. 


124.0 

-248.0 

126.4 


(A. 2) 


To  calculate  the  coefficients  f j ,  a  few  additional  geometrical  parameters 
(not  previously  defined  here  or  in  Figure  6)  are  required,  namely 


d  =14.33  cm 

w 

d  =15.04  cm 

c 

A^  =  11.218  cm^ 
I  ,  =  58.34  cm^ 

7.1 


(A.  3) 


The  parameter  S,  also  required,  is  defined  by  Equation  (33).  This 
parameter  thus  takes  on  the  value. 


S  = 


29,890  cm5 


(A. 4) 


The  coefficients  f^  are  also  functionally  dependent 
vertical  fixity  of  the  stiffener  ends.  Considering 
supported  case  leads  to  the  values, 


on  the  degree  of 
first  the  simply 


8.619 

-14.958 

6.958 


(A.  5) 
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The  coefficients  hj  and  f^  are  used  to  compute  the  three  coefficients  of 
the  quadratic  equation  in  K,  Equation  (49).  Using  the  above  values,  this 
results  in  the  expression 


165. IK  -  453. 3K  +  282.7  =  0 


(A.  6) 


for  which  two  roots  can  readily  be  computed,  namely 


Kl  *  1.788 
K2  =  0.958 


(A. 7) 


The  critical  tripping  loads  corresponding  to  these  two  roots  can  now  be 
readily  calculated  using  Equation  (46).  This  expression  requires  the 
vertical  moment  of  inertia  of  the  plate-stiffener  combination  which, 
through  the  routine  computation,  is  found  to  be 


I  =  2071  cm4 


(A. 8) 


The  functions  H^(K)  and  F^(K)  must  also  be  evaluated.  Using  the  two 
roots  for  K  and  the  coefficients  hj  and  fj  previously  defined,  the 
following  values  result. 


H1(K)  =  84.67  2.421 

F^K)  =  4.118  0.675 


(A.  9) 


Substitution  of  the  above  values  into  Equation  (46)  finally  leads  to  the 
desired  result,  namely 
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~  16, 540  N/cm 


for  K  =  K, 


2,885  N/cm 


(A. 10) 


for  K  =  K„ 


Of  these  two  values,  the  one  corresponding  to  the  root  K  Is 

“  ;;rst'  :he  crmcai  — *  a  of 

*  rr 

structure.  never  OCCur  in  the  real 

For  the  case  where  the  ends  of  the 
cedure  Is  the  sane  as  previous!,  descrlhed  ^ 

coefflcients  fj  t.ke  on  the  following  values, 


fQ  =  -0.381 
fj  =  3.042 
f2  =  -3.042 


(A. 11) 


These  values  and  the  copf f fr-iont-  . 

quadratic  equation,  ^  "l  ^  to  the 


and  the  roots 


369. 9K2  -  658. IK  +  282.7  =  0 


(A. 12) 


Kx  *  1.054 
K2  =  0.725 


(A. 13) 


These  roots.  In  turn,  lead  to  the  critical  tripping  lo,d. 
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f+  4,397  N/cm 


for  K.  =  K, 


(A. 14) 


-  38,040  N/cm 

for  K  =  K_ 

The  occurrence  here  of  tripping  loads  of  both  signs  (positive  and 
negative)  indicates  that  tripping  can  occur  due  to  loads  acting  in  either 
direction,  unlike  the  simply  supported  case.  Thus,  both  of  the  above 
tripping  loads  represent  realistic  tripping  phenomena.  Which  of  these 
loads  in  any  given  situation  is  the  one  of  practical  interest  will  depend 
on  the  direction  of  the  lateral  loading  or  loadings  under  consideration. 

FINITE  ROTATIONAL  RESTRAINT 

If  any  but  a  very  small  amount  of  rotational  restraint  along  the 
stiffener's  line  of  attachment  to  the  plating  is  to  be  considered,  then 
the  solutions  for  tripping  which  include  web  deformations  should  be 
utilized.  For  axial  end  loading  this  solution  is  given  by  Equation  (70), 
the  coefficients  for  which  are  provided  by  Equations  (71).  These 
expressions  involve  the  additional  parameters 


D  =  2.376  x  10  N-cm 

w 

z  =  11.26  cm 

2 

A  =  21.56  cm 
s 

I  =  880  cm4 

ps 


(A. 15) 


as  well  as  the  dimensionless  rotational  restraint  parameter,  R.  Assuming 
a  rotational  spring  constant  of  the  supporting  plating  of  C  =  30,000  N-cm/ 
rad  leads  to  a  value  for  R  of 


R  =  0.322 


(A. 16) 


For  this  value  of  R  the  coefficient  has  the  following  values  (for  m=l). 
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The  root  of  Interest  Is  the  second  one;  the  first  root  represents  a 
solution  mode  which  Is  mathematically  but  not  physically  possible. 

For  the  case  of  lateral  loading,  solutions  which  Include  the  effects 
ot  web  deformations  are  not  presently  available.  Consequently,  to  treat 
this  ease  one  has  the  choice  of  Inner  lug  the  Influence  of  the  rotational 
resistance,  or  It  this  course  Is  not  preferable,  of  making  use  of  an 
available  finite  element  program.  Finite  element  results  for  the  problem 
being  examined  here  are  presented  In  Table  5  In  the  main  t<'xt  for  several 
values  of  rotational  resistance. 

INKLASTIC  KITKOTS 

2 

If  the  material  has  a  yield  stress  ot  o  »  20,000  N/cm"  and  a 
structural  proportional  limit  ratio  p^  -  0.5  Is  assumed,  then  the  axial 
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end  loading  solutions  for  both  zero  and  finite  rotational  resistance  should 
be  adjusted  for  inelastic  effects  according  to  Equation  (100).  For 
Pr  -  0.5,  this  expression  becomes 


(<o  • 

ec 


re 


(°e)cre  1  2  °Y 


°Y 


4(o  ) 

e  ere  J 


(A. 20) 


V 


<°e>c 


re 


Thus  for  zero  rotational  restraint,  the  inelastic  tripping  stress  becomes 


(o  ) 

e  cr 


20,000 


20,000 

4  x  18,130 


*  14,480  N/cm 


2 


(A. 21) 


and  for  C  =  30,000  N-cm/rad 


(0  )  -  20,000 

e  cr 


20,000  ‘ 

4  x  34,700 


“  17,120  N/cnr 


(A. 22) 


It  the  "effective  width"  for  plating  effectiveness  under  inplane  loading 
Is  also  assumed,  in  this  example,  to  be  equal  to  30.5  cm  and  the  stiffener 
spacing  b  -  45  cm,  then  the  mean  inelastic  tripping  stresses  can  be 
computed  according  to  Equation  (27).  This  expression  gives 


(a  ) 
m  cr 


21. 5b  +  (30.5  x  0.8) 
21.5b  +  (45.0  x  0.8)  ' 


14,480 


“  11, 560  N/cm" 


(A. 2  1) 
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for  zero  rotational  restraint  and 


(o  ) 

m  cr 


21.56  +  (30.5  x  Q.8) 
21.56  +  (45.0  x  0.8) 


17,120 


*  13,670  N/cm 


2 


(A. 24) 


for  C  -  30,000  N-cm/rad. 

Since  the  classical  plate  buckling  stress  for  the  plating  in  this 
example  (a*160  cm,  b-45  cm,  t“0.80  cm)  according  to  Equation  (94)  is 

o  ,  -  o  .  =  7,880  N/cm2  (A. 25) 

pbe  pb 

it  should  be  immediately  clear  that  iterations  are  unnecessary  in  this 
case  and  that  the  appropriate  solution  (assuming  that  the  resistance 
represented  by  C  is  due  to  the  plating)  is  the  one  corresponding  to  C  =  0 
(zero  rotational  restraint). 

To  examine  how  the  iterative  procedure  is  carried  out,  consider  the 
same  stiffener  attached  to  plating  1.2  cm  thick  at  a  uniform  spacing  of 
38  cm.  In  this  case,  the  classical  plate  buckling  stress  (for  a=160  cm, 
b»38  cm,  t=1.2  cm)  is 

0  .  =  24,860  N/cm2  (A. 26) 

pbe 

Correcting  this  value  for  inelastic  effects  (using  pf  =  0.5  here,  also) 
gives 

o  ,  -  15,980  N/cm2  (A. 27) 

pb 

Since  the  dimensions  of  the  plating  do  not  enter  into  the  calculations 
for  (°e)cre  an<^  (og)cr,  t*1C  Prev*ous  values  computed  are  also  valid  here. 

I 


Since  C  and  (o  )  are  functionally  related  by  Equation  (95),  the  value  of 
6  c  r 

(oe)cr  calculated  for  C  -  30,000  N-cm/rad  can  be  used  as  a  "starting" 
value  in  the  iteration  process.  The  value  of  Cq  used  in  Equation  (95)  is 
defined  by  Equation  (96)  which,  in  this  case,  has  the  value 

Cq  -  158,100  N-cm/rad  (A. 28) 

2 

Substituting  the  critical  stress  (°e)cr  “  17,120  N/cm  into  Equation 
(95),  produces 

C  -  158,100  (l-  )  ■=  0 


or 


(A. 29) 


C  =  0 


This  value  may  now  be  used  as  the  second  estimate  of  C,  or  (hopefully)  to 
minimize  the  number  of  cycles,  a  value  may  be  selected  which  is  somewhere 
in  between  the  "assumed"  (0  30,000)  and  "calculated"  (OO)  values.  Try 
C  *  5,000  N-cm/rad  for  the  second  estimate.  This  gives 


(o  )  =  21,070  N/cm2 

e  ere 

(o  )  -  15,250  N/cm2 

e  cr 


(A. 30) 


and,  consequently, 

C  =  158,100  (l-  15?93~q)  =  7,220  N-cm/rad 
Now,  assume  C  ■  6,000  N-cm/rad.  This  results  in 


(A. 31) 
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(A. 32) 


(o  )  -  21,840  N/cm2 

e  ere  ’ 

(0  )  -  15,420  N/cm2 

e  cr 

C  =  5,540  N-cm/rad 


Now,  assume  C  =  5,800  N-cm/rad.  This  results  in 


(O  )  =  21,680  N/cm2 

e  ere 

(a  )  =  15,390  N/cm2 

e  cr 

C  =  5,840  N-cm/rad 


(A.  33) 


Because  of  the  very  slight  changes  in  the  tripping  stresses  (°e)cr 
between  the  last  two  cycles,  practical  considerations  would  dictate  that 
the  iteration  process  has  converged  sufficiently  for  engineering  purposes. 
Thus,  the  final  solutions,  rounded  off,  are 


(0  ) 

-  21,700  N/cm 

e  ere 

(0  ) 
e  cr 

=  15,400  N/cm' 

(A. 34) 


COMBINED  LOADING 

Returning  to  the  original  problem  (b=45  cm,  be=30.5  cm,  t“0.80  cm), 
consider  the  situation  when  both  axial  end  loads  and  lateral  loads  are 
occurring  simultaneously.  In  this  case  elastic  tripping  is  assumed  to 
occur  when  the  expression 


O 

e 


(o  ) 
e  ere 


+ 


q 

xre 


1 


(A. 35) 


is  satisfied.  If  rotational  restraint  is  ignored  (C-0) ,  then  values  of 
(og)cre  and  qcre  (for  clamped  end  conditions)  have  previously  been 
computed.  Thus,  the  above  expression  becomes 
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_ ®_  +  — 9 — 

18,100  A, 400 


1 


(A.  36) 


For  a  lateral  loading  of  magnitude  q  «  450  N/cm  (corresponding  to  a 

2 

uniform  pressure  of  10  N/cm  )  the  maximum  bending  stress  in  the  flange 

2 

of  the  stiffener  will  be  on  the  order  of  5,000  N/cm.  Thus,  at  tripping, 

Og  will  be  the  predominant  loading.  Solving  for  CTg  and  designating  the 

result  by  (o  )  to  distinguish  it  from  (o  )  ,  gives 

J  e  cre,q  e  ere 


(a  )  =  18,100  (l- 

e  cre.q  V  4,400  /  ^ 

-  16,250  N/cm2 

Since  a is  the  predominant  loading,  the  inelastic  correction  is  applied 

to  (o  )  resulting  in 

e  cre,q  ® 


(a  )  -  13,850  N/cm2  (A. 38) 

e  cr,q 


The  mean  axial  tripping  stress  is  then 

(a  )  -  11,060  N/cm2 

m  cr,q 


(A. 39) 
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APPENDIX  B 

QUADRATIC  COEFFICIENTS  FOR  H  (K)  AND  F  (K)  FUNCTIONS 

m  m 

The  functions  H  (K) ,  H  (K) ,  and  F  (K) ,  defined  in  the  main  text  in 
m  m  in 

relation  to  laterally  loaded  beams,  are  quadratic  functions  of  the  scalar 
parameter  K.  Thus  these  functions  can  be  written  as 

“.<»  *  h0  +  V  +  v2 

Hm(K)  =  hQ  +  h^K  +  h2K2  (B.l) 

Fm(K)  "  f0  +  fl(K)  +  f2K2 

The  determination  of  the  appropriate  value  of  K  requires  that  the  specific 
coefficients  h^,  h^,  and  f^  be  identified.  This  can  be  readily  accom¬ 
plished  simply  by  rearranging  and  collecting  terms  for  each  of  the  three 
functions.  Therefore,  from  Equation  (43a)  one  obtains 

2  2 

h.  =  (m+2) 2  +  (nri-2)4  (I  72+D  +  (B.2a) 

0  a2GJ  2  tt2GJ 

[2  2 
(m+2)2  +  (m+2)4  (I  72+D  +  (B.2b) 

a2GJ  2  tt  GJ  J 


h,  =*  2(m2+2m+2)  +  [m4+(m+2)4]  (I  s2+D  +  (B.2c) 

2  a2GJ  2  tt  GJ 

Similarly,  from  Equation  (43b),  the  coefficients  for  Hm(K)  are  defined 

h  -  (m+2)4  (^)  2  +  6 (m+2) 2  (1-v)  + 

0  V  a  ’  tt2D  d 


(B. 3a) 


-  2  f(m+2)4  (^)  +  6(m+2)2  (1_v)  +  3|fc_“J  (g 
L  tt2d  d 

W  __ 


2  (m  +(m+2)  ]  (~)  +  12(m2+2m+2)  (1-v)  +  ^ c 


7T2D  d 
w 


(B.3c) 


For  the  function  P.W.  two  sets  of  coefficients  must  he  defined,  one  for 

spy  supported  case  and  the  other  for  the  clamped  case.  For  the 
simply  supported  case  Equation  (45)  will  yieid 


f n  =  (®+2 ) 


0  o  bu  r 

£  _  -*  _ W  /  1  \  T 

it2  '  tt2s  L  ftf  (dc"h+  2  C)  ^  h 


fi  ”  ~  2(m+2)2  +  ■£-  l+m(m+2)  (l+  — l 

77  L  (  f«H. 


(m+1) ‘ 


12d  r 

+  7T  LAftf  (dc~h+  i c)  +bet2hj 

E2  “  ®2  +  (®+2)2  -  1+m (m+2)  (l+  — L_ 

77  L  l  (m+i  • 


12d 

'  77  IVf  (“s'*  I 0  +V2h] 

For  the  clamped  case  these  coefficients  become  (from 


(B.4a) 


(B.4b) 


(B.4c) 


Equation  (53)) 


•i  6d  r 

'°  ’  '  7  "  7i  [Vf  (de'b+  I  ')  +bet2|>] 


(B.5a) 
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12d 


tt2S 


Vf  (dc-h+  1  *)  +bet2h 


(B. 5c) 


(Note  that  in  the  expressions  for  the  f ^ ,  the  symbol  h  for  the  location  of 


the  neutral  axis  should  not  be  confused  with  the  coefficients  h^  previously 


defined. ) 
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APPENDIX  C 

COMPARISONS  WITH  EXPERIMENTAL  DATA 


Although  relevant  experimental  data  on  tripping  failure  are  virtually 
nonexistant,  the  Naval  Construction  Research  Establishment  (NCRE)  grillage 
tests  previously  cited  provide  data  on  two  steel  grillages  whose  ultimate 
failure  was  ascertained  to  have  been  caused  primarily  by  longitudinal  stif¬ 
fener  tripping.  The  pertinent  scantlings  (mean  values)  of  the  two  gril¬ 
lages,  identified  as  grillages  la  and  lb,  were  as  follows: 


Grillage 

la 

a  = 

48  in. 

b  = 

24  in. 

t  = 

0.315 

in. 

d  = 

6.05 

in. 

t  = 
w 

0.284 

in. 

f  = 
w 

3.11 

in. 

- 

0.56 

in.* 

*Nominal  values 


lb 

48  in. 

24  in. 
0.310  in. 
6.0  in.* 
0.28  in . * 

3.0  in.* 

0.56  in.* 


Both  grillages  were  nominally  identical,  the  only  difference  being  the 
presence  of  a  15  psi  lateral  loading  on  grillage  lb. 


GRILLAGE  la 

Employing  a  modulus  E  =  30  *  10^  psi  and  a  Poisson's  ratio  V  =  0.3 
the  elastic  tripping  stress  for  grillage  la  according  to  Equation  (70) 
with  C  =  0  and  m  =  1  is 


(a  )  =  104.7  ksi  (C.l) 

e'cre 

Substituting  this  value  into  Equation  (100)  with  cr^  =  37  ksi  and  p^O.5 
leads  to  an  inelastic  axial  tripping  stress 

(Oe)cr  -  33.7  ksi  (C. 2) 
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Since  the  classical  elastic  plate  buckling  stress,  a  given  by  Equation 

(94)  is  about  18.7  ksi,  the  assumption  of  the  value  C  ■  0  is  thus  seen  to 

be  correct  and  no  iteration  with  regard  to  this  parameter  is  necessary. 

Therefore,  one  may  proceed  directly  to  the  calculation  of  the  mean  tripping 

stress  (o  ) 

m  cr 

The  calculation  of  the  mean  tripping  stress  requires  the  knowledge  of 
the  plating  effective  width  bg.  Since  the  grillage  at  tripping  has  been 
loaded  above  the  level  of  plate  buckling,  it  is  clear  that  the  plating 

3 

should  be  expected  to  be  somewhat  less  than  fully  effective.  Assuming 
the  following  expression 


b 

e 

b 


2 

8 


where,  in  this  case,  8  is  a  plate  slenderness  ratio 


(C.3) 


8 


s  IS 

t  v  E 


(C.4) 


leads  to  a  value  of  effective  width 

be  0.607b  =  14.6  in.  (C.5) 

Inserting  this  value  into  Equation  (27)  (with  (°e)cre  replaced  by  (Ofi)  ) 
in  turn  provides  the  desired  mean  inelastic  tripping  stress, 

(o  )  =  24.5  ksi  (C. 6) 

m  cr 

This  value  compares  quite  favorably  with  the  experimental  value  of  27.8  ksi 
measured  at  NCRE. 
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GRILLAGE  lb 

Grillage  lb  is  nominally  identical  to  la,  differing  only  through  the 
addition  of  a  15  psi  lateral  pressure.  For  this  grillage,  using  Equation 
(70)  with  C  =  0,  m  =  1,  and  the  scantlings  previously  listed,  the  elastic 
axial  tripping  stress  becomes 

(0  )  =  99.4  ksi  (C. 7) 

e  ere 

(Note  that  this  value  does  not  yet  include  the  effects  of  lateral  pressure 
and  differs  slightly  from  that  of  grillage  la  because  of  slight  differences 
in  the  tabulated  mean  values  of  the  various  scantlings.)  To  include  the 
effects  of  the  lateral  loading,  the  tripping  load  due  to  lateral  loading 
alone  is  computed.  Assuming  an  effective  breadth  of  50t  for  the  combined 
plate-stiffener  combination  and  then  substituting  the  appropriate  parameters 
and  coefficients  into  Equation  (46)  results  in  a  critical  lateral  loading 
for  clamped  ends,  C  =  0,  and  m  =  1, 

q  =  +  17,100  lb/in.  (C.8) 


The  applied  lateral  loading  of  15  psi  represents  an  effective  line  loading 
of  360  lb/in.  for  24  in.  frame  spacing.  Using  the  linear  interaction  Equa¬ 
tion  (84)  the  elastic  axial  tripping  stress  in  the  presence  of  lateral 
loading  can  now  be  computed, 

(0  )  =  (o  )  (l  -  -3—  ^ 

e  cre.q  e  ere  ^  qcre  ) 

~  59-4  (l  -  ITTIoo)  <c'9> 

=97.3  ksi 
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Since  it  is  easily  demonstrated  that  at  failure  the  axial  load  is  predomi¬ 
nant,  the  inelastic  tripping  stress  can  now  be  calculated  in  the  same 
manner  as  for  grillage  la,  using  Equation  (100).  This  gives 

(a  )  =  33.5  ksi  (C.10) 

e  cr,q 

Assuming  the  same  expression  for  b^/b  as  for  grillage  la,  leads  to  a  value 
of  effective  width 

b  ss  0.600b  =  14.4  in.  (C.ll) 

e 

and  finally  to  the  mean  inelastic  axial  tripping  stress 

(a  )  =  24.1  ksi  (C.12) 

m  cr ,  q 

The  corresponding  experimental  value  from  the  NCRE  tests  is  27.1  ksi. 

Again  agreement  appears  quite  reasonable. 

Although  the  agreement  between  analytical  and  experimental  results 
appear  quite  encouraging,  it  is  only  fair  to  point  out  that  several  param¬ 
eters  whose  values  are  difficult  to  define  precisely  can  influence  the 
predicted  analytical  values  to  varying  degrees.  The  most  significant, 
perhaps,  is  the  choice  of  plating  effective  width  which  is  necessary  to 
determine  "mean"  from  "peak"  stresses.  It  can  be  demonstrated  that,  de¬ 
pending  on  the  degree  of  plating  effectiveness  assumed,  a  rather  wide  range 
of  mean  tripping  stress  values  can  result.  For  example,  in  the  case  of 
grillage  la,  the  inelastic  axial  tripping  stress  (ue)cr  =  33.7  ksi,  which 
is  a  peak  stress,  can  produce  mean  tripping  stress  values  (°m)cr  as  high 
as  33.7  ksi  for  fully  effective  plating  to  as  low  as  10.2  ksi  for  complete¬ 
ly  ineffective  plating.  The  most  appropriate  value  lies  somewhere  in  be¬ 
tween  these  extremes.  Indicating  that  such  wide  ranges  of  predicted  values 
can  exist  should  not  be  interpreted  as  completely  negating  the  value  of  the 
excellent  agreement  described  here,  rather  it  serves  to  point  out  that 
caution  must  be  exercised  in  attempts  at  experimental  correlation  when 
limited  amounts  of  data  are  available. 
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